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Abstract. 

A many-body theory of paramagnetic Kondo insulators is described, focusing 
specifically on single-particle dynamics, scattering rates, d.c. transport and optical 
conductivities. This is achieved by development of a non-perturbative local 
moment approach to the symmetric periodic Anderson model within the framework 
of dynamical mean-field theory. Our natural focus is the strong coupling, 
Kondo lattice regime; in particular the resultant 'universal' scaling behaviour 
in terms of the single, exponentially small low-energy scale characteristic of the 
problem. Dynamics/transport on all relevant (oj,T) scales are considered, from the 
gapped/activated behaviour characteristic of the low-temperature insulator through 
to explicit connection to single-impurity physics at high to and/or T; and for optical 
conductivities emphasis is given to the nature of the optical gap, the temperature scale 
responsible for its destruction, and the consequent clear distinction between indirect 
and direct gap scales. Using scaling, explicit comparison is also made to experimental 
results for d.c. transport and optical conductivites of Ce^BiiPtz, SmB 6 and YbBi 2 . 
Good agreement is found, even quantitatively; and a mutually consistent picture of 
transport and optics results. 

PACS numbers: 71.27.+a Strongly correlated electron systems; heavy fermions - 
75.20.Hr Local moment in compounds and alloys; Kondo effect, valence fluctuations, 
heavy fermions 
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1. Introduction. 

In the field of strongly correlated electrons, lanthanide- or actinide-based heavy electron 
materials constitute a longstanding challenge to experimentalists and theorists alike 
[1,2]. The majority of such systems, heavy fermions, are of course metallic; whether 
they be paramagnetic or ordered, Fermi- or non-Fermi liquids. Among them however 
resides a class of materials with insulating ground states: the so-called Kondo insulators, 
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containing a large variety of compounds reviewed e.g. in [3-7], and well known examples 
including Ce^Bi^Ptj,^ SmB 6 , YbBi 2 , CeRhAs, FeSi. This mainly cubic class of 
paramagnetic systems exhibit narrow-gap insulating/semiconducting behaviour at low 
temperatures, while their 'high' temperature behaviour is largely indistinguishable 
from metallic heavy fermions, amounting in essence to a lattice of / ions that scatter 
conduction electrons independently via the Kondo effect. The insulating gap has long 
been argued (see e.g. [3]) to arise from hybridization between essentially localized /- 
levels and a broad conduction band, the essential physics involving a flat /-band crossing 
one conduction band such that there are exactly two electrons per unit cell ('half-filling'); 
albeit that the resultant hybridization gap is not of course a simple one-electron entity, 
being strongly renormalized by many-body interactions that reflect the localized and 
hence correlated nature of the /-levels. As such, the Kondo insulators provide [3] a 
realisation of the simplest, canonical model for understanding heavy electron systems 
[1,2]: the half-filled periodic Anderson model (PAM), in which each lattice site contains a 
non-degenerate, correlated /-level hybridizing locally with a non-interacting conduction 
band, and which represents the natural lattice generalization of the single-impurity 
Anderson model (AIM) [2]. 

In the present paper we consider the half-filled, symmetric PAM within the powerful 
framework of dynamical mean-field theory (DMFT, reviewed in [8-11]). Formally 
exact in the large- dimensional limit, DMFT provides a tangible approximation in finite 
dimensions, whereby electron dynamics become essentially local but remain wholly non- 
trivial [8-11]. Our basic aims here are to provide a many-body description of dynamical 
and transport properties of paramagnetic Kondo insulators, specifically single-particle 
dynamics, dynamical conductivities and static electrical transport; and to develop the 
theory to the point where quantitative comparison with experiment can be made. 

These goals are of course easier stated than achieved, and the PAM has been studied 
extensively within DMFT via a wide range of techniques. Numerical methods include the 
numerical renormalization group (NRG) [12,13], quantum Monte Carlo (QMC) [14-17] 
and exact diagonalization [18]; while theoretical approaches include perturbation theory 
in the interaction strength [19,20], iterated perturbation theory [21,22], the lattice non- 
crossing approximation [23,24] and the simpler average t-matrix approximation [25], 
large- N mean-field theory [26,27], and the Gutzwiller approach [28,29]. NRG aside 
however, the above techniques suffer in general from well recognised limitations; whether 
it be an inability to handle large interactions and hence recover the exponentially small 
scales that are the hallmark of strongly correlated systems, failure to recover Fermi 
liquid behaviour at low-energies, unrealistic confinement to the lowest energies, and so 
on. Within DMFT all correlated lattice-fermion models reduce to an effective quantum 
impurity hybridizing self-consistently with the surrounding fermionic bath [8-11], i.e. 
to an effective, self-consistent AIM. Techniques for the latter thus underpin the former. 
Motivated in part by this we have been developing a 'local moment approach' (LMA) to 
quantum impurity models (AIMs) [30-35], the main emphasis of which is on dynamics 
and transport. Intrinsically non-perturbative and able to capture the spin-fluctuation 
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physics characteristic of the strongly correlated Kondo regime, the LMA encompasses 
all interaction strengths U and recovers simple perturbative behaviour in weak coupling 
[30]. Dynamics on all energy scales are handled, and the low-energy dictates of Fermi 
liquid behaviour satisfied (although the approach is not confined to Fermi liquid ground 
states, see e.g. [33,34]). Results for dynamics arising from the LMA have been shown 
[30-32,34] to give very good agreement with NRG calculations; and, for static magnetic 
properties, with exact results from the Bethe ansatz [35]. 

In a recent paper [36] we have further developed the local moment approach to 
encompass T = single-particle dynamics/spectra of the symmetric PAM. While plain 
perturbative behaviour is again recovered in weak coupling, the natural focus of [36] 
was on the strong coupling {i.e. large-C/) Kondo lattice regime. At sufficiently low 
energies in the vicinity of the Fermi level, the LMA recovers correctly the 'insulating 
Fermi liquid' behaviour [36] that reflects adiabatic continuity in U to the non-interacting 
limit of the simple hybridization-gap insulator [2,10]. This is manifest in preservation 
of the single-particle gap, now characterised by a renormalized gap scale A g which 
is reduced from its non-interacting counterpart by the quasiparticle weight Z that 
embodies many-body interactions. In agreement with consensus [12-15,24,28,29], strong 
coupling dynamics were found [36] to be characterized by the single low-energy gap 
scale A g ; which is exponentially small in strong coupling [36] (reflecting its dependence 
on Z), thus leading to a clear separation between low- and high-energy scales. In 
consequence, 'low'-u; dynamics exhibit scaling: being dependent solely and universally 
upon oo/Ag, with no dependence on the 'bare' high-energy parameters {U etc) that enter 
the PAM Hamiltonian. The simplest manifestation of scaling is that at sufficiently low 
energies oj/ A g the spectral behaviour amounts [2,10] to a quasiparticle renormalization 
of the non-interacting hybridization-gap insulator; which is of course the justification 
for renormalized band structure ideas. By itself however such quasiparticle behaviour 
gives rather a crude caricature of the scaling spectra, for it is confined to the immediate 
vicinity of the Fermi level uj = [36]; beyond which and on scales on the order of A g 
itself, non-trivial dynamics rapidly sets in, embodied in long, slowly varying spectral tails 
that reflect genuine many-body scattering/lifetime effects. As we show in the present 
paper, it is in fact this that dominates both dynamics and transport properties for all 
temperatures. 

The paper is organised as follows. The model and basic underlying theory is 
discussed in §2; formulated for an arbitrary host lattice, and including relevant aspects 
of the LMA in general (§2.1) as well as the specific class of diagrams contributing to the 
associated dynamical self-energies S ct (c<j;T) that we employ here in practice. §3 deals 
briefly with some basic formal results for conduction electron scattering rates, electrical 
transport and optical conductivities. Here and thereafter we consider explicitly and 
together two canonical host lattices [8-11], the hypercubic and Bethe lattices; our aim 
throughout being to emphasise both the differences and, more importantly, similarities 
between these two representative cases. Results arising from the LMA are then presented 
systematically in §s 4ff. Our primary emphasis is again the strong coupling Kondo lattice 
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regime of the PAM, this being both where the theoretical difficulties lie and the regime 
generally applicable to small-gap Kondo insulators. By the same token we focus largely, 
albeit not exclusively, on the scaling behaviour of dynamics and/or transport, now 
depending universally on uj = uj/A 9 andT = T/A g . This is important for many reasons, 
not least because the lack of scale separation inherent to some previous approaches 
has we believe led to a number of misconceptions in the literature, particularly in 
regard to the scales relevant to the T- and a>-dependence of the conductivity. In §'s 
4,5, and considering all Cj and T scales, single-particle dynamics, conduction electron 
scattering rates and d.c. transport are considered; including explicit connection to single 
impurity physics at high frequencies and/or temperatures ('Kondo logs' etc). Optical 
conductivities cr(co>; T) are considered in §6, with particular emphasis given to the nature 
of the optical gap and the clear separation between indirect and direct gap scales in both 
the ^-dependence and thermal evolution of cr(u;T). 

In §7 we turn to experiment, considering three prototypical materials for which 
extensive and reliable data is available [3-7], viz Ce^Bi^Pt^, SmB e and YbBi 2 ; our aim 
being direct comparison between the present theory and experimental results for both 
d.c. transport and optical conductivities. That may be achieved in a minimalist fashion, 
employing directly the scaling behaviour discussed in §s 4-6; which requires neither 
multiparameter fits nor in general a specification of the bare model parameters. Good 
agreement between theory and experiment is found, even quantitatively; with many of 
the characteristic features arising theoretically being directly apparent in experiment, 
and a mutually consistent description of transport and optics thereby arising. The paper 
concludes with a brief summary. 

2. Model and theory. 

In standard notation, the Hamiltonian for the PAM is given by: 

(i,j),(7 i,a i,a 

The first term describes the uncorrelated conduction(c) band with nearest neighbour 
hopping tij = t, rescaled as t oc t*/y/Z c in the large dimensional limit where the 
coordination number Z c — > oo [8-11] (with t* the basic unit of energy). The second 
term refers to the /-levels with site energies 6f and on-site repulsion U; such that 
€f = — y for the particle-hole (p-h) symmetric PAM considered here, for which 

nf = E CT < flfia >= 1 and n c = E CT < c t c ia >= 1 ( for a11 u ) as appropriate 
to the Kondo insulating state. The final term in equation ()2.1j) describes c//-level 
hybridization via the local matrix element V, whence the model is characterized by two 
independent dimensionless parameters, vizU '/i* and V/t*. 

Our natural focus is on local single-particle dynamics, embodied in the retarded 
Green functions G u (u)) («-»■ —i0(t) < {fia(t), fj a } >) and likewise G^(u)) for the c-levels, 
with corresponding local spectra D^(u) = --ImG^w) (and v — c or /). A knowledge 
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of local dynamics and their thermal evolution is in turn sufficient within DMFT [8-11] 
to determine optical and transport properties, as detailed in JJJ 

We begin with some brief remarks on the trivial limit V = 0, where (equation (|2.1jl ) 
the /-levels decouple from the free conduction band. The latter is specified by its local 
density of states po( e ) = _ -Ini<7o( e )> an d it proves useful in the following to denote by 
H(z) the Hilbert transform 



OC 



H(*) = / de EM (2.2) 

J — oo Z e 

for arbitrary complex z. The free c-electron (local) propagator go(uj) is itself given 
simply by 

</oM = tt{cu + ) (2.3a) 
= [u+ - Soiu)]' 1 (2.3b) 

with u + = uj + i0 + here and throughout; where (as used below) equation (|2.36|) defines 
the Feenberg self-energy Sq(uj) [37,38], with Sq(ui) = S[go] a functional of go (since 
go = H(S + l/g )). While our subsequent discussion holds for an arbitrary conduction 
band po(e), specific results will be given in §'s Aff for the Bethe lattice (BL) and 
hypercubic lattice (HCL); for which within DMFT the Po( e ) are respectively a semi- 
ellipse and an unbounded Gaussian, given explicitly by [8-11] : 

Po(e) = ^r[l-(eA*) 2 ^ : |e| < U BL (2.4a) 

A>(e) = —7— exp (-[e/t,] 2 ) HCL (2.46) 

As noted in ^ we are interested in the homogeneous paramagnetic phase of the 
PAM, for which the G^(uj) = G v (uo) are site- independent. The major simplifying 
feature of DMFT is that the /-electron self-energy is site-diagonal [8-11], and from 
straightforward application of Feenberg renormalized perturbation theory [37,38] the 
G v (uS) are given by: 

-1 

(2.5a) 



G f ( 



OJ 



V z 

' , V' 
uo — Tif (u; T) 



2 1 -1 



(2.56) 



u+-S(lo)_ 

Here £/(a;;T) is the conventional single self-energy (defined to exclude the trivial 
Hartree contribution which identically cancels ej = — j), such that £/(u;;T) = 
(u\ T) — i Sj(u;; T) with ^(u; T) > for all (uj, T); and with p-h symmetry reflected 

in 

S / (^T) = -[S / (-^;T)]* (2.6) 

together with S(lu) = - [S(-u)\* , G v (uj) = - ^(-o;)]* and hence D"(u) = D v (-uj) 
for the spectra. Equation ()2.56|) embodies the connection to an effective impurity model 
that is inherent to DMFT [8-11], for it may be cast in the 'single impurity' form G*(&) = 
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[uj + — T) — A cff (c<j)]^ 1 ; with an effective hybridization A eS (u) = V 2 [u + — S(u)] 1 

which, in contrast to that for a pure Anderson impurity model (AIM) and by virtue of 
its dependence on the Feenberg self-energy S(u), depends implicitly on coupling to the 
other sites in the correlated lattice and as such must thus be self-consistently determined. 
Specifically, the Feenberg self-energy S(u>) = S[G C ] is precisely the same functional of 
G c (uj) as it is of go(uj) in the V = limit (e.g. S = \t 2 G c for the BL). In consequence, 
G c is given using equation (J2.5aj) by 

G c (uj) = H( 7 ) (2.7) 



where 



V 2 

7 H =lu + - — — — — (2.i 



uj+ - S/(cu;T) 

(and we add in passing that in physical terms 7i(u;) = Im7(w) gives the conduction 
electron scattering rate which will be considered further in JJJ). 

It is this that, for an arbitrary po(^), prescribes the conventional 'single self-energy' 
route to the propagators G u (uj): given £/(a;; T), and hence 7(0;), G c (u) = B.(j) follows 
directly by Hilbert transformation; S(ui) follows (from equation (|2.5ajl ) as 

SM = 7-gfrt (2 - 9) 

and G*(u) then follows in turn from equation 1)2.56)1 . In practice of course the 
problem must be solved iteratively and self-consistently, because the approximate 
Tif(u;T) employed will itself in general be a functional of self-consistently determined 
propagators. Self-consistent second order perturbation theory in U [20] provides a 
direct example; as too does iterated perturbation theory (IPT) [21] where Sj(cj;T) 
is constructed from the second order (in U) skeleton diagram, employing host/medium 
/-electron propagators Q(oj) = [(G^(c<;)) _1 + S/(cu; T)] given from equation ()2.56)1 by 

-1 



g(u) 



v 2 



(2.10) 



lu+-S(lu)_ 

and thus dependent on S(u) itself. Results arising from IPT will be discussed in §'s 4-6. 



2.1. Local Moment Approach. 

In the conventional route to dynamics sketched above, the success of any particular 
theory is naturally determined by the approximation employed for the single 
self-energy £/. Therein lie well known difficulties [2], notably the inability of 
conventional perturbation theory to handle strong interactions in general, and to recover 
exponentially small scales that are the hallmark of strongly correlated behaviour; 
together with the divergences that plague perturbation theory in U [2] if one attempts 
to perform essentially standard diagrammatic resummations (e.g. RPA) of the sort 
one intuitively expects should be required to capture the regime of strong electron 
correlations. For these reasons the LMA [30-36] eschews an approach based directly 
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on the single £/, and focuses instead on a two-self-energy description that is a natural 
consequence of the mean-field approach from which it starts. 

There are three essential elements to the LMA. (i) First that /-electron local 
moments ('/i'), viewed as the initial effect of interactions, are introduced explicitly 
from the outset. The starting point is thus broken symmetry static mean-field (MF, 
i.e unrestricted Hartree-Fock); containing two degenerate, local symmetry broken MF 
states, corresponding to /i = ±|/i|. Grossly deficient by itself (see [36]), MF nevertheless 
provides a starting point for a non-perturbative many-body approach, to which end 
(ii) the LMA employs the two-self-energy description that is a natural consequence 
of the underlying two local saddle points. The associated self-energies are built 
diagrammatically from the underlying MF propagators, and include in particular a 
non-perturbative class of diagrams ( §2.21 and figure H below) that capture the spin- 
flip dynamics essential to describe the strongly correlated regime, (iii) The final key 
idea behind the LMA at T = is that of symmetry restoration [30-36]: self-consistent 
restoration of the broken symmetry endemic at pure MF level, and hence recovery on the 
lowest energy scales of the Fermi liquid/quasiparticle behaviour that reflects adiabatic 
continuity in U to the non-interacting limit. 

As detailed in [36] the G u (u), which are as they must be rotationally invariant, are 
expressed formally as (cf equations (2.5)) 

G» = |[qM + G[H] (2.11) 



where 



V 2 

uo + = S(u) 

lo + -E a (u;T) 

V 2 

u + - T, a (u; T) 



(2.12a) 
(2.126) 



LO+ - S(io)_ 

(and cr = t/| or +/— ); and with the /-electron self-energies separated as 

t a (u; T) = -|£7|/Z| + E a (oj; T) . (2.13) 

The first term here represents the purely static Fock bubble diagram that alone 
is retained at pure MF level (with given explicitly by equation ()2.17|) below). 
The second term E CT (o;;T) = T,^(u;T) — i£^(u;;T), is the all important dynamical 
contribution mentioned above, with p-h symmetry reflected in 

E (T (a;;T) = -[S_ ff (-u;;T)]* (2.14) 

(such that G v a (u) = - [G v _ a {-uj)]* and hence G v {u) = - [G"(-u)}*). 

Equations (|2.11H 2) are the two-self-energy counterparts of the single self-energy 
equations (2.5). And for an arbitrary conduction band Po(e) and given {S^u;; T)}, they 
may likewise be solved straightforwardly (cf the above discussion of equations 1)2.71 -9)): 
defining 

V 2 . . 
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such that G c [uj) = | [7^ — S] 1 (equations (|2. 111 12a)). and comparing to G c {uj) = 
[7 — S}~~~ 1 (equations (|2.5ol 8)). the 7 CT 's are related to the single 7(0;) (equation (|2.8jl ) 
by 

[i(7 T H-7iH)] 2 



7(w) = 3 [7tM +7iM] + 



j[7 T (w) + 7iH]' 



(2.16) 



Given £ CT (u;;T) and hence 7t(w), this equation together with S{u) = 7 — 1/H(7) 
(equation (|2.9jl ) may be solved iteratively for S(u) (employing an initial 'startup' S, 
say S = jtlgo(u)); which procedure is both straightforward and numerically fast. And 
with S(uj) then known the G v (u) follow directly from equations (|2. 11( 12). In particular, 
the underlying MF propagators may be obtained from this procedure in one shot, the 
static MF self-energies being given by S CT = S° = —ax with x = \U\p,\. For any given 
x, the MF propagators g^{oj) and hence spectra d v a {uj) = d^(u;x) thus follow; and at 
pure MF level the local moment \fi\ is then found from solution of = \p,\, where the 
Fock bubble |/Z| = \ft(x)\ appearing generally in equation (|2.1Hjl is given by 



du 



dUu)] f{u-T) 



(2.17) 



with f(uj;T) = [e W//T + l] the Fermi function. 

The single self-energy Hf(uj;T) likewise follows as a byproduct of the above 
procedure, since solution of equations (|2.9|) . (|2.16|) given {S cr (c<j;T)} determines both 
S(lj) and 7(0;), whence (see equation ()2.8|) ) Sy(u;;T) = uo + — V 2 [u + — 7(u;)] _1 thus 
follows; which relation may be recast equivalently as 



+ 





~\ ^E|(o;; T) — E;(o;; T)j 




2 




S T (w;T)+E;(w;T) 



(2.18) 



where Q{oS) = oo + — V 2 {uo + — S{uo)) 1 is precisely the host/medium /-propagator 
equation (|2.10|) . The resultant conventional single self-energy may thus be obtained 
directly, given the two self-energies S CT (o;;T) equation (|2.1H|I (although not of course 
vice versa); and the particular class of diagrams contributing to the dynamical S cr (c(j; T) 
that we retain in practice are specified in §2.21 

As mentioned above and discussed further in [36], the final, central idea behind 
the T = LMA is self consistent restoration of the broken symmetry inherent at MF 
level. This is embodied mathematically in the symmetry restoration (SR) condition 



E T (u; = 0;T = 0) = Sj 



UJ 



0; T = 0) at the Fermi level to — 0; and hence E a (0; 0) = 



(for either a) for the p-h symmetric PAM here considered, 



i.e. 



Sf(0;0) = Ef(0;0) - \U\p\ = 



(2.19) 

(using S CT (0;0) = S^(0;0)). In physical terms, satisfaction of SR ensures [36] 
that single-particle dynamics on the lowest energy scales amount to a quasiparticle 
renormalization of the non-interacting limit U — 0, reflecting Fermi liquid behaviour 



Dynamics and transport properties of Kondo Insulators 



9 



in the general sense of adiabatic continuity to that limit. For U = (= £/) the non- 
interacting Green functions are denoted by g$(uj; V 2 ), with spectra d^ui] V 2 ) and the 
^-dependence explicit; such that (via equations (|2.2[ 3.5.7.8)) d^u; V 2 ) = p (u> — V 2 / u) 
and d Q (u\ V 2 ) = ^d^u; V 2 ) with po(e) the free (V = 0) conduction band density of 
states, e.g. equations (2.4). For U = and all V ^ the system is thus of course a simple 
hybridization gap insulator [39], with a gap A° g (V 2 ) that is hard for the generic case of a 
bounded Po(e) (e.g. the BL, equation (|2.4ajl ) and (strictly) soft for an unbounded po(e) 
satisfying po(e) — > as |e| — > oo (e.g. the Gaussian characteristic of the HCL, equation 
(12.46)0 . If the SR condition equation (j2.19|) is satisfied then (see [36]) for U > the 
leading, lowest-u; behaviour of the full T = G v (uj) follows as 

G c (u) ~ g c (u; ZV 2 ) (2.20a) 
G f (u) ~ Zgl(uj;ZV 2 ) (2.206) 

where Z = [1 — (dT, R (u>; Q)/du)) U} =o]~ 1 is the quasiparticle weight (given equivalently by 
[36] Z = [1 — (d'E R (uj;0)/du) ul =o}~ 1 ); resulting in preservation of the insulating gap, 
A g = A®(ZV 2 ), that is generically reduced from the non-interacting hybridization gap 
by the quasiparticle weight factor Z, with Z <C 1 in strong coupling. 

Equations (2.20) embody the quasiparticle behaviour of the PAM on the lowest 
energy scales, i.e. the 'insulating Fermi liquid' behaviour that evolves continuously from 
the non-interacting limit. Imposition of SR equation (|2.19|) . as a single condition at the 
Fermi level u = 0, naturally underlies the LMA; and amounts in practice, as detailed 
in [36], to a self-consistent determination of the local moment \p\ (superseding the pure 
MF condition \p\ = \p(x)\, see equation (j2.17)0 . 

2.2. Dynamical self-energies. 

Our final task is to specify the class of diagrams retained in practice for the dynamical 
T)'s (equation (2.13)). These embody self-consistent dynamical coupling of single- 
particle excitations to low-energy transverse spin fluctuations and are precisely as 
considered in [36] for T = 0, extended to finite-T following arguments identical to 
[32] for the Anderson impurity model; the reader is thus referred to [32,36], for full 
details. The diagrams are summarized in figure ^ where wavy lines denote the local 
interaction U, the double line propagator denotes the broken symmetry host /medium 
/-electron propagator Q^ a (oS) specified below (equation ()2.21|l ). and the local /-level 
transverse spin polarization propagator is shown as hatched; diagrammatic expansion 
of the resultant in terms of MF propagators and dynamical self-energy insertions is 
discussed in [36,40]. 

In physical terms the diagrams shown in figure Q describe correlated spin-flip 
scattering processes that are essential to capture in particular the strong coupling Kondo 
lattice regime of the model: in which having say, added, a a-spin electron to a — a-spin 
occupied /-level on site i, the — cx-spin hops off the /-level thus generating an on-site 
spin-flip ( embodied in the transverse spin polarization propagator); the — cr-spin electron 
then propagates through the lattice/host in a correlated fashion, interacting fully with 
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Figure 1. Class of diagrams for the /-electron self-energies here retained 

in practice. The interaction U is denoted by a wavy line and the renormalized 
host/medium propagator (see text) by a double line; the transverse spin polarization 
propagator is shown hatched. 



/-electrons on any site j ^ i (reflected in the host /medium G before returning 
to site i at a later time whereupon the originally added a-spin is removed (and which 
process simultaneously restores the spin-flip on site i). 

The renormalized medium /-propagator G- a (uj), which embodies correlated 
propagation of the — a-spin electron through the lattice, is given explicitly by [36] (c/ 
equation (EHU^ : 

V 2 



2.21 



LU+ - S(u) 

As used below, Q a {uS) may be separated as Q a {uj) = G^iyj) + G^{w) into the one-sided 
retarded Hilbert transforms, given by 

G a H = / dux — — — (2.22) 

J-oo UJ~UJ l + l0+ 

with T> a {uj) = — -1mQ(r{ijj) the corresponding spectral density (and 6(x) the unit step 
function). 

Specifically, the retarded LMA £j(u;;T)(= — uj; T)]* by p-h symmetry) is 
given by 

S T (c;T) = f/ 2 f"** [°° dw 2 X + -("i; T) V ^ hfa U2 ) (2.23a) 

with 

h(u i; u 2 ) = 0(wi)/(a; 2 ; T) + 0{-u x ) [1 - f(u 2 ; T)} (2.236) 

(which reflects the hard core boson character of the local /-level spin flips [32], 
whose statistics are dictated by the probability with which fermions can hop from/to 
site i to/from the surrounding host lattice, as embodied in the Fermi functions). 
X + ~(uj; T) > denotes the (local) spectral density of transverse spin excitations, given 
by x + ~(co>; T) = sgn(o;) Imn + ~(w; T) with U + ~(uj; T) the retarded, finite-T polarization 
propagator. The latter is given at the simplest level, considered here, by an RPA-like 
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p-h ladder sum in the transverse spin channel, obtained by straightforward analytical 
continuation of the imaginary time 

n+-(iu; m ) = °U + -(iu m ) [1 -U !^-^)]' 1 . (2.24) 

The bare polarization bubble diagram °II +_ may itself be expressed either in terms of 
the broken symmetry MF propagators {gf(uj;x)} ( fl2.1j) . as shown explicitly in figure 
lc of [36] and with the resultant LMA referred to therein as LMA I; or in terms of 
the self-consistent medium propagators {Q a }, correspondingly referred to as LMA II. In 
practice as shown in [36], results for single-particle dynamics arising from LMA I /II are 
very similar, and for that reason explicit results are given in the present paper for LMA 
I alone. We also add that for T = 0, equations (2.23) reduce generally (via equation 

flZZfl ) to 

£ T (^; 0) = U 2 I — x + ~(^i; 0) kwO^Cw! + u) + + w)l (2.25) 

J — oo 7T L 

as employed (in time-ordered form) for T = in [36]. 

The above considerations specify the LMA two-self-energies that we consider in 
practice; X CT (co>; T) being given in its entirety by equation (|2.13jl . \p\ therein by equation 
(12.17)1 and the dynamical £-|-(u;; T) by equations (2.23) ( with T) by p-h symmetry). 
The problem is readily solved numerically. As explained in £ 12.11 (following equation 
()2.16Jl ). for given equations (|2.9U 6) may be solved straightforwardly for the 

Feenberg self-energy S(u) and the full c//-electron propagators G u (u). An iterative, 
self-consistent solution is naturally required, since the {S^} are functionals of the 
renormalized medium /-propagators {Qa}, themselves given explicitly by equation (|2.21j) 
and thus dependent on S(u). And since both Q a and the MF propagators g[ depend 
explicitly onx = ^U\[/,\, it is numerically optimal to solve for fixed x, with U determined, 
as opposed to vice versa [36]. 

The problem is first solved for T = 0, ensuring that symmetry restoration (equation 
(J2H2J)) is satisfied at each iterative step. For any given x = \U\fi\ solution of equation 
()2.19|) determines U, and the local moment = |/x(T = 0)| then follows directly. As 
explained in [36] this step generates a spin-flip resonance in x + ~(co>; T = 0) centred on a 
non-zero frequency u m . This is the low-energy scale characteristic of the Kondo lattice 
(and with uj m oc Z, the quasiparticle weight); its origins within the LMA thus stemming 
from SR, and its physical significance being that it sets the timescale r ~ h/uj m for 
restoration of the locally broken symmetry inherent at the crude level of pure MF. For 
T > the same iterative algorithm may be employed, except that SR is no longer 
required. Temperature enters the problem in two distinct ways: (a) explicitly, and 
centrally, via the Fermi functions inherent in h(ui, uj 2 ) and x + ~{^i] T) (equation (2.23)); 
and (b) implicitly, via the T-dependence of the local moment = \[i(T) \ in x = ^U\fj,\. 
The latter may be encompassed via |/i(T)| = |/i(0)| + S\fi(T)\, with |/x(0)| the T = 
moment required to satisfy SR as above; and with 8\fi(T)\ calculated in practice at MF 
level [32], such that 5|/x(0)| = 0. As for the Anderson impurity model [32], we find 
however that the resultant T-dependence of is negligible for essentially all (U,V) 
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(provided one is not concerned with physically irrelevant temperatures of the order of 
U), and we thus omit it from the results shown explicitly in §'s 4ff. 

3. Electrical transport and optical conductivity. 

Within the large-dimensional framework of DMFT a knowledge of single-particle 
dynamics, in particular the (uj, Independences of the /-electron self-energy Hf(u;T), 
enable q = transport properties to be determined [8-11]. In this section we specify 
some basic results for the conduction electron scattering rate and cj-dependent electrical 
conductivity, for both the hypercubic and Bethe lattices. These are independent of the 
particular approximation employed to determine £/(o;;T). But they naturally underlie 
the results obtained via the LMA that are given in §'s 5ff. 

Equation (|2.7J) for the conduction electron Green function, G c {uj) = H(7) with 
7(0;) = uj + — V 2 [uj + — S/(u;;T)]~ 1 , is equivalently but more familiarly expressed as 
G c {uj) = N' 1 J^al^ ~ e a — £ c (u;; T)] -1 . Here e a denote the states of the free (V = 0) 
conduction band with density of states (equation (2.4)) po(e) = N" 1 J2 a ^( e ~ e «) ( e -S- 
e « = £k fc> r a Bloch decomposible lattice); and S c (u;;T) (= uj + — 7(0;)) is the purely 
local conduction electron self-energy, related to the /-electron single self energy £/ by: 

S C (^;T) = V 2 [uo + - Hf^-T)]- 1 (3.1) 

It will prove useful in the following to rewrite equation (2.7) as 

/oo 
dep (e)G c (e;uj) = < G c (e; u) > e (3.2) 
-00 

with the e-resolved propagator G c (e;u) = [7(0;) — e] -1 = [u + — e — S c (u;;T)]~ 1 and 
corresponding spectrum D c (e; uj) = — ^ImG c (e; oj)\ and where <A(e) > t = f depo(e) A(e) 
defines the e-average of any A(e). 

In particular, the conduction electron scattering rate 1/t(u; T) {K = 1) considered 
in §4.1, is given by 

/ 1 ^ = 7j(^; T) = -ImS c (a;; T) (3.3) 
t(uj; 1 ) 

(with 7/ = Inry). It is conveniently expressed in the dimensionless form l/f(u;T) = 
7/(a;;T) = 7rpo7/(^; T), with p — Po( e = 0); and is given in terms of the /-electron 
self-energy £/ = — i by 

1 _ - r #n = ^M,T) 

r{^T) 7n ^' } [o/-A^ 1 Ef(a;;T)]" + [A^ 1 Ej( W ;T)]3 K ' } 

where uj' = 00/ A and A is defined by 

A = 7rV 2 p . (3.5) 

Physically, A is the hybridization strength that would arise for a pure Anderson 
impurity model (AIM), in which a single correlated /-level is locally coupled (via V) 
to the free metallic conduction band po{e). And equation (3.4) will prove important in 
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connecting results for the PAM at large uj and/or T, to those for the pure AIM itself 
(and notwithstanding the fact that the ground state for the latter is metallic, while that 
for the symmetric PAM is of course insulating); see also equation (3.12) below. 

Calculation of the cu-dependent conductivity is in principle straightforward in 
the large-dimensional limit of DMFT, since the strict absence of vertex corrections 
[41] to the (q = 0) current-current correlation function means that only the lowest 
order conductivity bubble diagram survives [8-11]. We denote the trace of the 
resultant conductivity tensor by a{uj;T) (| of which, denoted by a(u;T), provides an 
approximation to the isotropic conductivity for a d = 3 dimensional system). This may 
be cast in the form 

M = F„( Wi r) (3.6) 

where o"o = ] ^jf-y — 7^ such that a ~ 10 4 — 10 5 [Ocm] -1 for lattice constants a in 
the physically realistic regime of 1 — 10 A. The dimensionless F a (u; T) depends on the 
lattice type, and is given explicitly for the hypercubic lattice (HCL) and Bethe lattice 
(BL) by 

,2 roc 

F hcl (uj; T) — — dw x [fM - f{u x + u)\ < D c (e; u^D^e; u 1 + u) > e (3.7a) 

V J -00 

t 2 f°° 

F bl {uj;T) =^ dui [fM - f(ui + u)} < D c (e; ui) > e < D c (e; ui + u) > e (3.7b) 
v J -00 

(with f(uj) here the Fermi function). The exact result equation (3.7a) for the HCL is 
of course well known (see e.g. [8-11]) and widely used even in studies employing the BL 
(e.g. [21]); although for the latter we emphasize that it is equation (3.7b) which follows 
from direct analysis of the conductivity bubble diagram. Equations (3.7a,b) correspond 
in an obvious physical sense to limiting forms of behaviour, from fully correlated to 
uncorrelated averages of the D c (e; u;)'s. Both will be considered in §'s 5ff. 

Before proceeding we comment briefly on evaluation of Fhcl(w,T) itself, 
under a single approximation: namely that in < D c (e; ui)D c (e; cu 1 + uj) > t = 
J^lo depo(e) D c (e; ui)D c (e; uj\ + uj) entering equation (3.7a) the e dependence of po(e) 
is neglected, Po( e ) — Po(0) = Po- A straightforward integration over e e (—00,00) then 
yields 

7T [7a(wi +u)- 7M^i)] + [7/(^i + + 7/(^i)J 
(where 7(0;) = 7_r(w) + i7/(w)). Use of this in equation (3.7a), relating 7 to the 
conduction electron self energy S c as above, gives 

F HCL (u;T)~ Im / dw x — ■ — r— — — 3.9 

■kuj J_ OCJ uj - S Cjr (^i + uj) + £ c , a (u;i) 

where r/a here denote retarded/advanced functions; and for the particular case of the 
d.c limit uj = 0, equation (3.8) with equations (3.4,7a) give: 

Fhcl(0]T)~ l -[pvQ 2 J duj ' 9 ^ f(uj;T) = l -[p,Q 2 <f> (3.10) 
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Equations (3.9,10) are likewise well known [25] and widely used; and the latter in 
particular, expressing the d.c conductivity in terms of an averaged scattering time, 
is physically intuitive. We emphasize nevertheless that they are approximate (granted 
even the legitimate neglect of vertex corrections), being dependent on the 'flat-band' 
approximation po(e) ~ PoVe as above; and it is in fact simple to show that this 
approximation by itself fails to account for the existence of the Kondo insulating gap 
that is characteristic of the symmetric PAM. That said however, one expects physically 
that equation (3.10) should provide a good approximation to the d.c conductivity at 
least for sufficiently high temperatures dominated by incoherent scattering; the question 
here, considered briefly in §5, being how high is high? 

Finally, for explicit use in §5, we consider the pure Anderson impurity model itself, 
denoting by Pi mp (T) the change in resistivity due to addition of the impurity to the non- 
interacting host, and P; mp (T) = Pim P (^ 1 )/pimp(0). This has the same functional form as 
equation (3.10), viz (see e.g. 



1 /"^ Jp^ fimp(a;; T) (3.11) 



P'imp( T ) J-oo 9uJ 

The impurity scattering rate is given by [2] l/fi mp (o;;T) = 7rAo-Di mp (u;; T) with 
D imp (uj;T) the impurity spectral function and A the hybridization strength equation 
(3.5) (such that 7rA Am P (0; 0) = 1 from the Friedel sum rule [2]). Denoting the impurity 
single self-energy by £(u;; T) = T, R (u; T)— i £/(c<j; T), the impurity transport rate is given 
explicitly by 

f imp (^;T) [u'-A^(u;T)}^+[l + A^(u;T)}^ l ' ' 

with uj' = u/A ; which should be compared to its counterpart for the PAM, equation 
(3.4). 



4. Single-particle dynamics. 

We turn now to single-particle spectra resulting from the LMA specified in §2. Our 
focus for obvious physical reasons is the strong coupling (large-C/) Kondo lattice regime, 
wherein universality arises; in particular the resultant scaling behaviour of dynamics, 
and the thermal destruction of the Kondo (insulating) gap. Conduction electron 
scattering rates l/f(u;T) are also considered here (§4.1) since these are closely related 
to single-particle dynamics (as directly evident from equation (3.4)). 

In strong coupling, as shown in [36] for T = and consistent with previous work 
[12-15,24,28,29], the problem is characterised by a single low-energy scale; as embodied 
in A g given by 

V 2 

A g = Z— (4.1) 

with Z — [1 — (i9Sj(cj; 0)/du) U! =o] ~ 1 the quasiparticle weight. This sets the scale for 
the Kondo gap, which (§2.1) is generically hard (as for the Bethe lattice), albeit strictly 
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soft for the hypercubic lattice; other embodiments of this low-energy scale, such as 
the spin-flip u m arising in x + ~( UJ ] 0) (§2.2 and [36]) are equivalent to A g , being simply 
proportional to it. There are two distinct issues relating to the low-energy scale. First 
its dependence on the 'bare' high-energy parameters entering the Hamiltonian, viz the f- 
electron Coulomb repulsion U, hybridization V and bandwidth scale t* (or equivalently 
p 1 oc t*). In strong coupling it is known, from NRG results in particular [13], that the 
gap scale becomes exponentially small; such behaviour is indeed found within the LMA 
(viz lnA g oc — U/(8V 2 p ) as detailed in [36]). This is important insofar as it guarantees 
a clean separation between low-energy and/or temperature scales on the order of the 
Kondo gap and multiples of it; and high u and/or T on the order of the bare energy scales 
(Hubbard satellites in the /-electron spectra to cite an extreme). Failure to recover this 
pristine separation of scales can obscure much relevant physics in the strong coupling 
regime (we provide examples in the following sections). This may arise either because 
strong coupling and/or low temperatures are difficult to access in a numerical approach 
(e.g. QMC); or because the approximate theory used leads to an insufficiently small 
low-energy scale, for example algebraic rather than exponential decay in U . Iterated 
perturbation theory (IPT) [21,22,10] provides an example of the latter, and is discussed 
further in the following sections. 

Granted a clean separation of energy scales however, the precise dependence of 
the gap scale A g on bare parameters is subsidiary. Of primary importance is that 
physical properties, whether dynamics or transport, exhibit universal scaling behaviour 
on experimentally relevant 'low' uo and/or T scales on the order of the Kondo gap 
and (in principle arbitrary) multiples thereof. The scaling is of course in terms of 
the gap A g itself, and is thus independent of the bare model parameters; in contrast 
to the correspondingly non-universal behaviour arising on truly high w or T scales 
characteristic of the bare parameters. Neither is such scaling of purely theoretical 
interest, for an experimental gap of (say) 50K ~ 4meV is tiny on the eV scale typical 
of bare parameters, and many multiples of it may arise before non-universal scales are 
reached in practice. Moreover it is arguably less preferable, as well as unnecessary in 
general, to have comparison to the tv and/or T dependence of experiment hinge on a 
particular choice of bare model parameters; as opposed to a knowledge solely of the 
experimental gap scale, which alone is required if comparison is made instead to the 
relevant scaling form (§7). 

Figure 2 summarizes salient results for the T = /-electron spectrum [36]; showing 
irA Df(u) (with A = irV 2 p , equation (3.5)) as a function of Co = uo/ ' A g [42], and for 
both the BL (solid line) and HCL (dashed). This is the universal scaling form arising 
from the LMA [42], with no dependence whatever on any of the bare parameters U, 
V, t*. It naturally depends on the lattice type, but only weakly and on scales up 
to a few times the gap, beyond which the two scaling spectra coincide. The Kondo 
insulating gap is directly apparent in figure 2. For the BL we also show explicitly 
the limiting quasiparticle form, whose recovery at sufficiently low-d; embodies adiabatic 
continuity to the non-interacting limit ('insulating Fermi liquid' behaviour, §2.1 and 
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Figure 2. T = scaling spectra nA Df(uj) vs Co = ui/A g for the BL (solid line) and 
the HCL (dashed line). The limiting low energy quasiparticle form is also shown (for 
the BL, dotted line). 



[36]); given explicitly using equation (2.20b) by irA D f (oj) ~ (4/|u>| 2 ) [1 - 1/uj 2 ]* 
for > 1 (and zero for |u)| = |u;|/A ff < 1, the gap). For \Co\ > 3 however the 
quasiparticle form is simply inadequate: it decays rapidly as ~ 1/|^| 2 and fails to 
recover the long, slowly decaying tails evident in figure 2. The latter, which dominate 
the scaling spectra at (moderate to) large Co — and in consequence transport properties 
at correspondingly 'high' temperatures (§5) — are logarithmically slow and discussed 
further in the following sections. Universal scaling in terms of Co is not of course confined 
to the /-electron spectra: the T = c-electron spectrum t*D c (uo) (or equivalently 
D c (uj)/pq) behaves likewise; and in consequence, as follows straightforwardly using 
equations (2.5), Aq 1 S/(cj;0) also scales universally (which is why the dimensionless 
conduction electron scattering rate has been defined as in equation (3.4)). 

For finite temperatures, what one expects in scaling terms is obvious; viz that the 
f/c spectra and /-electron self-energy should now exhibit universal scaling in terms 
of both Co = co/Ag and T = T/A g . That such scaling correctly arises within the 
LMA is demonstrated in figure 3 where, for fixed T — 2, the /- and c-electron scaling 
spectra are shown for three different interaction strengths U = U/t* = 5.6,6.1,7.0 and 
V 2 /tl — 0.2. The inset shows n AqD^ (co) on an 'absolute' scale, i.e. vs oo/t*, illustrating 
the exponential reduction in the gap scale with increasing interaction strength. The 
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main figures by contrast show the //c-spectra vs oj = cu/A g , from which the scaling 
collapse is evident (and in practice sets in by rather moderate interactions U ~ 4 or so). 
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Figure 3. Universal scaling spectra (BL) for fixed T = T/A g = 2 : n Ao D f (u) and 
D c (u)/p vs u>/A g for V 2 /tl = 0.2 and U = U/U = 5.6 (dashed), 6.1 (dotted), and 
7.0 (solid). Inset: corresponding /-spectra on an absolute scale, vs u/t*. 



The scaling illustrated above arises generically, and figure 4 shows the resultant 
LMA scaling spectra (for the BL) for a range of scaled temperatures up to T = 25. The 
thermal destruction/infill of the Kondo insulating gap is evident, occurring as expected 
physically for temperatures on the order of the gap A g , and accompanied in the case 
of the /-electron spectra by a corresponding destruction of the characteristic T = 
spectral 'horns'. We also note, as evident from the inset to the /-spectrum which shows 
the spectrum on an enlarged Co scale, that for any T the high-frequency asymptotics of 
the scaling spectra coincide with that for T = 0; as likewise expected physically, and 
which behaviour arises for frequencies |u)| ^> T. 

To give an 'all scales' perspective on thermal evolution, figure 5 shows the /- and 
c-electron spectra for the BL (specifically for U = 6.1 and V 2 /tl = 0.2) on an absolute 
energy scale, vs u>/t*; and for temperatures T = T/A g = 0, 10, 50, 450. The /-spectrum 
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Figure 4. Temperature dependence of the BL scaling spectra: 7rA -D^(u>) and 
D c (oj)/p n vs uj/Ag for temperatures T = T/A g = (solid), 1 (short dash), 2 (double 
point-dash), 5 (dotted), 10 (long-dash) and 25 (point-dash). Inset: /-spectra on an 
enlarged scale, out to Co — 100. 



is naturally dominated by the non-universal Hubbard satellites at ~ ^ which carry 
almost all the spectral weight (and are of course 'projected out' of the scaling spectra); 
the key low-energy universal spectral features shown in figures 2-4 being nigh on invisible 
in figure 5 as expected, since their net spectral weight is of order Z <C 1. The c-spectrum, 
which shows weakly remnant Hubbard peaks, is by contrast dominated by the envelope 
of the free conduction band (the |u;| < t* semiellipse for the BL); and again the thermal 
destruction of the low-energy spectral gap is barely visible. We add moreover that until 
temperature reaches essentially non-universal scales, the T-dependence of single-particle 
dynamics is confined to the relevant low-energies illustrated in figures 3,4 (see e.g. the 
inset to figure 5 where the Hubbard satellites are enlarged, noting that T = 450 here 
corresponds to T ~ 0.2i*). 
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Figure 5. Spectra on all scales (BL with U/U=6.1 and V 2 jt\ = 0.2): nA D f (u) and 
D c (uj)/p n vsuj/t* for temperatures T = T/A g = (solid), 10 (dashed), 50 (point-dash) 
and 450 (dotted). Inset: Hubbard satellites on a much reduced scale. 



4-1- Scattering rates. 

The conduction electron scattering rates T _1 (a;;T) are now considered. These are 
directly related to the /-electron self-energy Aq 1 Ej(cj;T) as in equations (3.3,4); and 
determine the conductivity via the dependence thereof (equation (3.7)) on the c-electron 
spectrum D c (e;uj) (or rather more directly via equation (3.10) in the latter's regime of 
applicability, considered explicitly in §5). 

Before proceeding we note that our conclusions regarding the scattering rates 
differ significantly from the work of [21]. In order to explain typical experimental 
conductivities for Kondo insulators, it was argued in [21] that scattering rates in the 
vicinity of the Fermi level should be on the order ~ 0.1 — 1 of the bandwidth t* (values 
some 2-3 orders of magnitude higher than for a clean metal like Cu). Scattering rates 
calculated in [21] were however found to be some 1-2 orders of magnitude lower than 
required. While the authors of [21] note that this behaviour is very surprising, they 
attribute it to an intrinsic limitation of the model itself; rather than e.g. to a limitation of 
the approximate calculations employed, or to a possible misidentification of the relevant 
temperature scales involved. We find by contrast, as shown below, that scattering rates 
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can certainly attain values on the order of the bandwidth t*; and indeed argue that rates 
of this order must arise over a wide, experimentally relevant temperature regime. 



0.8 



^* 0.6 

P 
3 

' P 0.4 



0.2 



0.5 



0.25 



-100 




100 




co = co/A 



g 



Figure 6. Thermal evolution of the conduction-electron scattering rates: r 1 (ui;T) 



(in units of the BL bandwidth t* 



w/A 3 for temperatures T — T/A g = 0.1 



(solid), 0.5 (dotted), 1 (short dash), 2.5 (point dash), 5 (long dash) and 20 (double 
point-dash). Inset: on an expanded scale, out to Cj = 100. 



The uj- and T-dependences of the scattering rate arising in strong coupling are 
illustrated in figure 6: viz r _1 (u;; T) in units of the bandwidth t*, v s Co = u/ A g and for 
a range of temperatures T = T/ A g from 0.1 to 20; as given explicitly via equation (3.4) 
(where the 'bare' factor of uJ = Cu^ plays no role since oc Z is exponentially small in 
strong coupling). The inset to figure 6 shows the same results on an expanded scale out 
to Co = 100 (showing that for \Co\ ^> T the 'tail' behaviour reduces to that characteristic 
of T = 0, the logarithmically slow decay of which is considered explicitly below, figure 
8). The only remaining ^-dependence arises on non-universal (and essentially irrelevant) 
scales on the order of the Hubbard satellites \u\ ~ where as expected physically the 
scattering rate is strongly peaked [21]; this is omitted from figure 6 for obvious reasons. 

The first point to note about r~ 1 (c<j;T) is that for T = precisely it contains a 
(^-function contribution at the Fermi level u> = 0; specifically nt*5(Co) as follows generally 
via equation (3.4) using Sj (u; 0) ~ — ^ as uo — > 0. For T ^ this becomes the resonance 
evident in figure 6, which naturally broadens with increasing temperature; and the only 
thermal scale on which this can occur is the gap A g — the single low-energy scale 
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characteristic of the system in strong coupling. This argument is rather general. The 
scattering rate in the vicinity of the Fermi level, |u)| < T, in consequence diminishes with 
increasing T from essentially arbitrarily large values (reflecting the insulating nature of 
the T = state); and does so on temperature scales of the order of the gap, A g . 
The results of figure 6 show moreover that for temperatures T = T/A g in the range 
~ 1 — 20, scattering rates in the relevant Cj regime are indeed on the order o/~ 0.1 — 1 
of the bandwidth t*. We also add that while a temperature range of this order certainly 
encompasses that relevant to experiment (§7), the essential point is not dependent on it, 
since for T 3> 1 we find the Fermi level scattering rate to decay slowly with temperature 
(specifically t _1 (cj = 0;T) oc l//n 2 (T), see equation (5.7) below). 

The issue of scattering rates was considered in [21] using IPT [43], explicitly so for 
a particular choice of bare parameters U/t* — 3 and V/t* = 0.25; and at a temperature 
T = O.lt*, a significant fraction of the free conduction bandwidth and some 50 times 
the corresponding IPT gap A g = ZV 2 /t :¥ . The resultant scattering rate in the vicinity 
of the Fermi level was found to be r _1 (0;T) ~ 10~ 2 t*, with which the authors of [21] 
support their view mentioned above. We have also performed IPT calculations for the 
same U and V; and indeed for T = O.lt* recover the results of [21]. We have further 
investigated the T-dependence of IPT over a wide T = T/ A g range (as well as a broad 
(U, V) range). Significantly, we find that for temperatures T up to ~ 5, the behaviour of 
the resultant t~ 1 (lu\T) is qualitatively similar to that shown in figure 6: in particular, 
scattering rates in the vicinity of the Fermi level are again found to lie in the range 
O.lt* — t*. For temperatures T > 3 — 5, the IPT scattering rates decay much more 
rapidly with increasing T than those arising from the LMA. This reflects the inability of 
IPT to capture the logarithmically slow decays in u> and/or T that are characteristic of 
the model in the strong coupling/scaling regime, as illustrated below (figure 9); and is in 
turn related to the fact that IPT leads to an algebraically rather than an exponentially 
small gap scale in strong coupling, and hence does not produce a 'clean' separation 
of low/high energy scales. Nonetheless if it was indeed the case that non-universal 
temperature scales on the order of e.g. O.lt* were pertinent in relation to experiment, 
then the resultant scattering rates would in general be too small to explain observed 
conductivities. Our view is that this is not the case, but rather that the relevant thermal 
scale for comparison to experiment is the gap A g and multiples of it; and on which scales 
the transport rates readily attain values on the order of ~ O.lt* — t*. Further support 
for this view will be provided in the following sections. 

We turn now to more detailed consideration of the conduction electron scattering 
rate in strong coupling, in particular its high-frequency 'tail' behaviour evident in figure 
6 (inset) and its relation to the /-electron single-particle spectrum. We begin with the 
latter. It is straightforward to show generally (using equations (2.5) together with the 
definition equation (2.8) of 7(0;) and equation (3.4)), that the asymptotic behaviour of 
the dimensionless scattering rate f _1 (o;;T) (equation (3.4)) in fact coincides with the 
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/-spectrum, viz 



(4.2) 



This holds asymptotically for f~ l {u]T) <C 1; which in practice means (|u;|/A ff =) 
| cZ> | ^> 1 for any T (the 'tails') or T > max(l, \oj\) for any \uj\. This is illustrated in 
figure 7(a) for T = 0, where the cu-dependence of n AqD^ (oj) and f _1 (u;; T) are compared 
(omitting the S(u) contribution for clarity): they naturally differ very significantly at 
low frequencies, but for |u)| > 10 or so their tails rapidly become coincident. The same 
holds at finite T as illustrated in figure 7(b) for T = 2.5 and 10; and for T 3> 1 the two 
coincide asymptotically for all oj. 




Figure 7. ttAqD* (oj) (solid lines) and scattering rates f 1 (ui;T) (dotted lines) vs 
uj/A g (for the BL). Figure 7(a) is for f = T/A g = 0, and figure 7(b) for f = 2.5 and 
10 (in an obvious sequence). 



An obvious question arises: what is the functional form of the |d)| ^> 1 'tails' 
in f~ 1 (cj;T) (or equivalently ttAqD^uj)), which as seen in figure 6 (inset) arise at 
sufficiently high frequencies for any temperature T? For large |cD| and/or T one expects 
[1,2] on physical grounds that properties of the Kondo lattice such as the /-electron 
spectrum or the resistivity p(T) (§5), should asymptotically approach those of an 
Anderson impurity model (AIM) (i.e. the pure Kondo model in strong coupling). The 
high-frequency behaviour of the impurity single-particle scaling spectrum D imp (ui) for 
the AIM has only recently been uncovered, using the LMA [31] (which gives excellent 
agreement with NRG results [34] for that problem). On the testable assumption that 
the high-frequency behaviour of the /-electron self-energy for the PAM has the same 
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Figure 8. T = scattering rates f 0) vs cj/A g for the BL (solid line) and HCL 
(dashed line) . Explicit comparison is also made to the asymptotic behaviour equation 
(4.4) (dotted line). 

scaling form as that for the AIM, we thus expect Aq 1 Sj(c<j) to behave as [31] 

Ao'SiH ~ ^[1 + ^ln 2 (a\Co\)] (4.3a) 

J 3 7T 

16 

A^Ef (u) ~ -sgn(w)—ln(a\u>\) (4.36) 

for \oj\ ^> 1 (or ^> max(l,T) at finite-T), with a a pure constant 0(1); and hence from 
equations (3.4) and (4.2) that 



1 37T 2 


[/n 2 (a 






f(u;T) 16 [/n 2 (a 
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+ [7r/n(a 


^l)] 2 



(4.4) 



In figure 8 for T = 0, the cu-dependence of f _1 (a;;0) is shown, for both the BL and 
HCL. It is also compared directly to equation (4.4). The asymptotics for the two lattices 
coincide in practice for |u)| > 5 — 10 (as expected e.g. from figure 2). And the 'high'- 
energy behaviour of the tails in t _1 (cj; T) (or equivalently tiAqD^ (a;)) is indeed seen to 
be that of equation (4.4) (with the constant a ~ 0.55 determined numerically). We also 
add in passing that the full form equation (4.4) is required for the very good agreement 
over a wide cu-range evident in figure 8, i.e. it is not exclusively dominated by the 
ultimate large-u) asymptotic behaviour nA D-f(u) ~ f _1 (u;;T) ~ 37r 2 /[16/n 2 (|u>|)]. 
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Figure 9. Comparing the w/A 9 -dependence of scattering rates and /-electron spectra 
(BL) arising from the LMA and IPT: 7rA £^(w) (LMA: dashed line; IPT: double 
point-dash line) and f _1 (u;;0) (LMA: solid line; IPT: dotted line). The LMA results 
are naturally universal; IPT results were obtained for U/t. t = 3 and V/t* = 0.25. 

Finally, explicit comparison to IPT [21,10] is made in figure 9, for the (Independence 
of both the T = /-electron spectrum 7rA -D^(u;) and the scattering rate f _1 (u;;0) 
(again omitting the 5(0) contribution). The LMA results shown are naturally the 
universal scaling forms, independent of the bare model parameters. Corresponding 
IPT results are again obtained for U/t* = 3 and V/t* = 0.25. The figure illustrates 
three points, (i) That IPT does not capture the important logarithmic tails (and hence 
e.g. produces much reduced scattering rates, even for T = considered), (ii) Its /- 
electron spectrum recovers, but amounts to little more than, the limiting low-frequency 
quasiparticle form (as seen by comparison to figure 2 above), (iii) The approach does 
not produce an exponentially small gap scale, and hence a 'clean' separation between 
low (meaning universal) and high energy scales; e.g. non-universal scales on the order of 
uj ~ O.lt* are reached by Cj ~ 50 in the IPT results shown in figure 9. We also note that 
the IPT transport rate r _1 (c<j; 0) actually vanishes for \uj\ < 3A 9 , rather than A g itself. 
This in turn is related to the spurious behaviour in the T-dependence of the resistivity 
arising from IPT, that is discussed further in the following section. 
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5. D.C. transport. 

We now consider d.c. transport properties of the PAM. Our aim here is to understand the 
temperature dependence of the static conductivity/resistivity, in particular over the full 
T-range in the strong coupling/Kondo lattice scaling regime; as well as the connection 
between the high-T behaviour of the PAM and that of the single-impurity Anderson 
model. 

As discussed in §3 the d.c. conductivity cx(0;T) = |<7o-F a (0; T), with the 
dimensionless F a (u;T) given respectively for the hypercubic and Bethe lattices by 
equations (3.7a,b). In the strong coupling regime of interest the static conductivity 
naturally exhibits scaling in T = T/A g , with no dependence on the bare system 
parameters U, V, t*. The resultant scaling form does of course depend on the lattice 
type (cf §4 for single-particle dynamics), and F a (0;T) for the two lattices is shown 
in figure 10. For the HCL the figure also shows (dotted) the T-dependence arising 
from the approximate but physically intuitive form equation (3.10), viz Fhcl(0;T) oc 
~~ f-tt ^^TkT^O^' T) with f(u; T) the (dimensionless) scattering time. Two initial 
points should be noted here. First that the latter approximation, while qualitatively 
inadequate as T — > where it fails to capture the insulating ground state, is very good 
for T > 2 or so; in fact as further evident from the inset to figure 10 where comparison is 
made up to T = 200, it coincides asymptotically, and rapidly, with Fhcl{0] T) obtained 
from the exact equation (3.7a). Second, recalling (§3) that the constant a is realistically 
on the order of ~ 10 4 — 10 5 (ficm) _1 , we note that the results of figure 10 readily 
encompass static conductivities a(0;T) on the order of ~ 10 3 — 10 4 (f2cm) _1 that are 
typical of Kondo insulating materials around room temperature [3-7] (see also §7). 

The insulating nature of the ground state is self-evident in figure 10; and for 
sufficiently low T one expects and finds activated transport of form 



which defines the transport gap A tr . For the BL we find A tr = A g , while for the HCL 
one naturally finds A tr oc A g , but with a proportionality constant differing from unity 



We turn now to consider the T-dependence of the static conductivity/resistivity over 
the full T-range, for both the HCL and BL; and in addition to compare the resultant 
behaviour, at high-T in particular, to that arising for the AIM. The resistivity for the 
pure impurity model, p' imp (T) = p imp (T) / p imp (0) with p imp (T) the change in resistivity 
due to addition of the impurity, is given from equation (3.11). To compare like with 
like in the following, comparison of equation (3.10) for the HCL with equation (3.11) 
for the AIM suggests that for the HCL we consider the dimensionless resistivity 



<t(0;T) oc exp[-A tr /T] 



(5.1) 



(~ 0.4). 




_ \[poQ 2 



(5.2) 
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Figure 10. Temperature dependence of the static conductivities in the Kondo lattice 
scaling regime: F a (ui = 0; T) vs T — T/A g for the HCL (solid line) and BL (point-dash 
line); the approximate equation (3.10) for the HCL is also shown (dotted line). The 
high-temperature asymptote for the BL (4/7T 2 ~ 0.4) is indicated by an arrow. Inset: 
results on an expanded scale out to T — 200. 



First let us consider the important differences between the hypercubic and 
Bethe lattices; which in turn reveals equally important similarities between the two. 
From figure 10 it appears that with increasing T, F HC l(0',T) is growing apparently 
unboundedly, albeit relatively slowly; while Frl(0; T) is slowly saturating to a constant 
value Frl(0; T) 0.4. This is indeed the case, and stems from the fact that in the limit 
of vanishing hybridization V = — where the /-levels decouple from the conduction 
band — the BL has a non-zero residual (T = 0) conductivity/resistivity (in contrast of 
course to any V > 0, for which the T = conductivity vanishes); specifically, as follows 
straightforwardly using equation (3.7b), that: 



F^°(0;T = 0) = [p t^ = - 



IT' 



(5.3) 



It is this, the residual conductivity of the free (V = 0) conduction band, that corresponds 
to the ultimate high-T constant asymptote seen in the inset to figure 10 for the BL 
(4/-7T 2 = 0.405.., marked by an arrow in the Figure). This is further evident if one 
considers the temperature dependence of F BL (0;T) on 'all scales', i.e. vs T/t*; as 
shown in figure 11 for a sequence of increasing U/t* in strong coupling. Here the T /in- 
dependence of Fg^°(0; T) is also shown, and is seen to constitute an upper bound to the 
T/t*-dependence of Fbl(0;T) itself. Since the gap A g is exponentially small in strong 
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coupling, the scaling behaviour of F BL (0] T) illustrated in figure 10 for finite T = T/A g 
corresponds in figure 11 to exponentially small values of T/t* = f^ 3 -; whence as seen 
from figure 11, Fbl(0; T) is bounded above by F^°(0; T = 0) for large-T in the scaling 
regime (which corresponds formally to any finite T in the limit A g — > 0). 




0' 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 

T/t, 

Figure 11. Temperature dependence of static conductivities on 'all scales', vs T/t*: 
F BL (0;T) for V 2 /tl = 0.2 and £//i*=4.6 (solid), 5.1 (dotted), 6.1 (dashed) and 
7 (point-dash). The V = limit Fg^°(0;T) is also shown (double point-dash). 
Inset: Corresponding results for Fhcl(0',T) for U/t*=A.6 (solid), 5.1 (short dash), 
5.6 (dotted), 6.1 (long dash) and 6.6 (point dash); in this case, F^£(0;T) = oo. 

The latter behaviour is readily understood in physical terms, since in the scaling 
regime the gap in the conduction band D c (u) rapidly fills in with increasing T and 
approaches the free V = conduction band (see e.g. figures 4b, 5b). For T ^> 1 the 
gap A g is thus in essence irrelevant, i.e. might as well be zero — which is of course 
just the V = limit, with its characteristic residual resistivity for the BL. For the HCL 
by contrast the residual resistivity for V = is precisely zero; reflecting the infinite 
conductivity associated with the coherent Bloch states that in this case arise when 
V = 0. So for the HCL by the same reasoning, one does not therefore expect the 
conductivity in the scaling regime to saturate for T ^> 1, as indeed evident in figure 
10 (see also figure 11 (inset) where the T/t^-dependence of F HC l(0;T) is shown for a 
sequence of increasing U/t*, as well as figure 12 below). We believe this conclusion 
to be rather general: if for whatever reason the free (V = 0) conduction band is 
characterised by a non-vanishing residual conductivity/resistivity, be it due e.g. to the 
intrinsic nature of the free conduction band or to the presence of disorder, then we 
expect the conductivity in the scaling regime to approach asymptotically this limiting 
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value for T>1. 

To take the above into account when comparing the two lattices it is natural to 
subtract from the resistivity any finite high-T asymptote, as embodied in l/F a (0;T) — 
1/F^ =0 (0; 0) (a redundant operation for the HCL where F^q L {0; 0) = oo); considering 
therefore — in obvious parallel to the AIM — the 'change in resistivity due to coupling 
the /-levels to the host conduction band'. Specifically, for the BL we consider below 
the following dimensionless resistivity: 
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Figure 12. Resistivities in the Kondo lattice scaling regime: p'hcl^T) (equation 
(5.2), solid line) and p' BL (T) (equation (5.4), point-dash line) vs T = T/A g . Inset: 
same results on a logarithmic T scale. 

We have just highlighted the differences between the HCL and BL. But their 
important similarities are evident in figure 12, where we show the resultant scaling 
behaviour of the resistivities PhclC^) (equation (5.2)) and p' BL (T) (equation (5.4)) 
as a function of T — T/A g out to T ~ 300; the inset shows the same results on a 
log scale, indicating that the p'(T)'s ultimately vanish as T — > oo. Two immediate 
features are apparent. First that the high-T asymptotics of the p'(T)'s coincide 
(whence re figure 10, the rate at which the BL conductivity Tbl(0;T) asymptotically 
approaches its saturation value is the same as that with which its HCL counterpart 
grows unboundedly); the leading asymptotics is obtained explicitly below. In fact for 
T > 5 — 10 or so, the behaviour of the two lattices is barely distinguishable. Second, 
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cursory inspection of figure 12 might suggest that the resistivities p (T) are themselves 
plateauing with increasing T, since they change only slightly over a large temperature 
range out to many multiples of the gap. This interpretation, while rather natural at 
first sight, is not correct since the p (T) — > as T — > oo. But it is indicative of the very 
slow T-decays involved; and for which, see §7, we believe there is experimental evidence 
in Kondo insulating materials. 

We now obtain the leading high-T behaviour of the p (T), here considering explicitly 
Phcl(T). Precisely the same result (equation (5.8) below) is obtained for p' BL (T), which 
separate calculation is outlined in the Appendix. From equation (5.2) together with 
equation (3.10) (which is asymptotically valid at large- T as illustrated in figure 10), 

where the scattering time t(uj; T) is given by equation (3.4) in terms of the /-electron 
self-energy Aq 1 Ej(o;;T). The latter is given for \uj\ ^> max(l,T) by equation (4.3); the 
generalization of which to T ^> 1 and any |u)| (holding also for the BL) is [32] 



A^Ej(a;;r) ~ ^ln 2 [ay/\u\* + (vrf ) 2 ] (5.6) 

with A^Ef^T) oc [A x Ej(w; T)] 1 / 2 ( SU ch that it may be neglected in equation (3.4) 
for the scattering time). Equation (3.4) thus yields 

16 

f(u; T) ~ A^Ej( W ; T) ~ —ln 2 (T)L(y; T) (5.7) 

where y = u/f, and L(y; f) = [1 + ln\a^v 2 + y 2 ]/ln(f )f such that L{y;f) -> 1 as 
T — > oo for finite y. Using equation (5.7) in equation (5.5), and changing variables 
therein from uo to y, gives directly the leading high-T behaviour of p^ Ci (T) as 

p'(T)^ 1 -^ (5.8) 
^ V ; 16/n 2 (T) V ; 

(which as above also holds for p' BL (T), see Appendix). 

Equation (5.8) reflects the anticipated connection between the PAM at high-T 
and the Anderson impurity model; being also the exact high-temperature asymptote 
for the strong coupling AIM (with Kondo temperature oc A g ) [2], first obtained [44] 
for the Kondo model from the leading logarithmic sum of parquet diagrams. That 
connection is seen more generally in figure 13, where the full temperature dependence 
of the HCL resistivity is compared directly to corresponding LMA results for the AIM 
(analogous comparison for the BL is clear from figure 12). Specifically we show p^ CL (T) 
vs T = T/Ag (with A g = ZV 2 /t :t = A Z/y/ir); compared to the impurity resistivity 
p' imp (T) = Ptmp (T)/p imp (0) vsf = T/A g ™P (where = A Z imp /^ with Z imp the 

AIM quasiparticle weight). The LMA p imp (T) is detailed in [32] where it is shown 
to agree rather well with NRG calculations for the AIM [45], to be asymptotically 
exact at high-T (equation (5.8) being recovered); to agree well with the Hamann 
approximation [46] (obtained by further parquet resummation) in the latter's regime 
of applicability down to T ~ 1, and to cross over smoothly to the AIM Fermi liquid 
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form p' imp (T) — 1 oc — T 2 as f — > 0. As evident in figure 13 the PAM resistivity 
Phcl(T), which exhibits activated insulating behaviour (equation (5.1)) for T < 0.5 or 
so, progressively diminishes with increasing T and for T > 1 — 2 essentially coincides 
with that for the AIM. 
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Figure 13. Comparison of resistivities for the HCL and Anderson impurity model: 
p'hcl(T) (solid line) and Pi mp (T) (dashed line) vs T — T/A g on a log-scale. Inset: 
comparison of p' HCL (T) arising from the local moment approach (solid line) and IPT 
(point-dash line). 

The above results are in marked contrast to those arising from iterated perturbation 
theory [21,10], the behaviour of which is qualitatively wrong. IPT results for p HCL (T) 
vs T are compared to the LMA in the inset to figure 13. The IPT calculations were 
performed for U/t* — 3 and V/t* = 0.25 (as employed for the BL in §4.1); and the 
inability of IPT to recover an exponentially small gap scale means that non-universal 
temperatures on the order of T ~ O.lt* are reached by T ~ 60. It is evident from figure 
13 that the IPT resistivities are generally much less than their LMA counterparts; and 
no hint of logarithmic T-decays are evident (which mirrors the absence of logarithmic 
tails in the ^-dependence of the IPT single-particle spectra and scattering rates, §4). 
Most significantly however, we see that the IPT resistivity actually vanishes as T — > 0; 
rather than exhibiting the correct divergence symptomatic of the insulating ground 
state. At first sight this is surprising, since IPT is qualitatively correct in predicting 
a T = gap in the single-particle spectra D c ^(cu) [21,10]; we have thus examined the 
problem in some detail. Its physical origins stem from the fact that the imaginary part 
of the IPT T = self-energy, £i(u;; 0), is zero for \u\ < 3A ff rather than for |u;| < A g as 
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one might expect (where 'zero' strictly means exponentially small for the HCL with its 
soft-gap, although that is irrelevant in the following); and in consequence the scattering 
rate r" 1 (u;; 0) (see equation (3.4)) likewise vanishes for < \uo\ < 3A 9 . But for \uo\ > A g , 
there is a high density of conduction electron states embodied in D c (uo) (see e.g. figures 
2 or 4). Hence states in the interval A g < < 3A g are essentially unscattered; and 
it is this that, on initially increasing T from zero (and thus accessing states in this 
interval), leads to an arbitrarily large conductivity and hence vanishing of p' H cl(T) as 
T — > 0. This can be demonstrated numerically in several ways (e.g. by 'cutting off' 
the Fermi functions df(u;T)/ 'duo entering p' HCL {T) in the above |u;| interval); but the 
physical argument above is the basic origin of the problem. This spurious behaviour 
does not arise within the LMA, for which S^(cj;0) and the single-particle spectra are 
non- vanishing in the same \u\ intervals. The problem can be circumvented within IPT 
by ad hoc addition of a finite-77 in frequency factors u + = 00 + 17] (to mimic an additional 
white-noise inelastic scattering rate); as has been employed in a weak coupling IPT 
study [47]. In this case we find the ultimate low-T behaviour of the resistivity to be of 
form P HC l(T) ~ nexp{bA g /T) (with b ~ 0.4); but with a crossover temperature scale 
to such activated behaviour that is entirely determined by 77 and vanishes as r\ — > 0. 

6. Optical conductivities. 

The natural progression is now to consider the optical conductivities o~(uj; T) = 
^cr F a (uj; T)), with the F a (uj; T) given explicitly for the hypercubic and Bethe lattices by 
equations (3.7). Before proceeding to specific results arising from the LMA, we consider 
a basic question regarding the optical gap: what is it? While strictly an issue at T = 
since all gaps are technically destroyed for T > 0, this of course has major ramifications 
for both the frequency and temperature dependence of the optical conductivity. 

We note at the outset that our most basic conclusion here differs qualitatively from 
the work of [21,15]. It has been argued hitherto (see e.g. figure 2 of [21]) that the 
optical gap in a(uj; 0) corresponds to the minimum direct gap A^ r of the renormalized 
band structure, as opposed to the indirect gap A in d oc A g = ZV 2 /t* that is manifest 
in the single-particle spectra D c (u) or D^(u); a conclusion that in turn underpins the 
interpretation of experiment given in [21] (see also [6]). We do not believe this to be 
correct on general grounds, neither do we find it supported by the present theory. By 
contrast we show, in agreement with the qualitative conclusions of [7], that : (i) the 
optical gap corresponds to the indirect gap (albeit that the direct gap scale is obviously 
also manifest in the (^-dependence of the optical conductivity); and (ii) that it is the 
indirect gap which sets the temperature scale for the thermal evolution of the optical 
conductivity in the Kondo lattice regime, in particular the 'filling' of the optical gap with 
increasing temperature (which we regard as entirely natural, since in strong coupling 
A g is the characteristic low-energy scale of the system). 

It is first necessary to explain why we disagree in particular with the conclusions 
of [21]. As discussed in §2.1 (see also [2,10,36]) the limiting low-frequency behaviour of 
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the single-particle Green functions amounts to a renormalization of the non-interacting 
limit; arising from the exact equations (2.5) by neglecting the imaginary part of the 
/-electron self-energy Ej(cj;0) and replacing (_ReS/(c<j; 0) =) {uj; 0) by its leading 
low-cu behaviour, viz 

S / (cu;0)~Sf(^0)~-(|-lV (6.1) 

with Z the quasiparticle weight. This leads directly to the quasiparticle behaviour 
embodied in equations (2.20), e.g. G c (uo) ~ <?o(o>; ZV 2 ) with gl(u\ V 2 ) referring to the 
U — limit; whence the lowest-u; behaviour of D c (uo) (and likewise D- f (uj)) is a simple 
quasiparticle renormalization of the non-interacting limit, with V 2 — ► ZV 2 . In particular 
(§2.1) the gap in the single-particle spectra is generically preserved. This is the indirect 
gap scale defined by, 

V 2 

A ind = 2A g = 2Z— (6.2) 

which is indeed manifest in the single-particle spectra (the 2 simply reflecting the 'full' 
gap, as opposed to that measured from the Fermi level at mid-gap, uj — 0). The 
quasiparticle form for the single-particle propagators is not approximate: it must be 
satisfied asymptotically at sufficiently low frequencies, reflecting as it does Fermi liquid 
behaviour in the sense of adiabatic continuity to the non-interacting limit; and indeed it 
is, as illustrated in figure 2 (see also [36]). We emphasise however that this quasiparticle 
behaviour is confined to the lowest frequency scales \uj\ = |u;|/A fl , up to at most a few 
times the (indirect) gap A g oc Z\ beyond which, as detailed in §4 and [36], this simple 
picture no longer holds. 

The arguments above may be extended — but with a wholly different validity — to 
the optical conductivity; it is this that underlies the conclusion of [21]. From equations 
(3.7), the optical conductivities are determined by the D c (e; uj) = —^ImG c (e; uj); where 
(§3) G c (e;uj) = [7(0;) - e]" 1 with 7 (w) = cu + - V 2 [uj + - E/^T)] -1 . If the asymptotic 
form equation (6.1) for E/(o;;0) is again employed here, then ^(oo) xs uo + — and 
hence G c (e;cu) [lu + — — e] _1 which again amounts to a renormalization of 
the non-interacting limit. The two branches characteristic of this renormalized band 
structure then follow from the zeros of the approximate [G c (e;u)}~ 1 ; being given by 
oo±(e) — |[e ± \J e 1 + 4ZV 2 ] , with the e-dependent gap function A^(e) = ou + (e) — c<J_(e). 
The minimum such gap occurs for e = 0; this is the (approximate) direct gap, 

A dir ~ 2^ZV (6.3) 

(which we add in passing is radically different from the indirect gap equation (6.2) in 
strong coupling where Z < 1, since A dir /A ind oc Z^ 1 ^ 2 ). If now one interprets the 
optical conductivity of the hypercubic lattice (given from equation (3.7a)) in terms of a 
naive picture of renormalized interband transitions, then since A dir above corresponds 
to the minimum energy for such, one would clearly expect it to vanish for uj < A dir ; 
concluding [21] therefore that the optical gap coincides with the direct gap, while the 
indirect gap is by contrast manifest only in the single-particle spectra, see e.g. figure 2 of 
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[21]. And indeed, if the above approximation G c (e; uj) ks [uj + — ^ — e] 1 is used explicitly 
in equation (3.7a), one finds directly that Fhcl(w; 0) = for all uj < A dir = 2^[ZV 
(and in fact likewise for Fhcl(w, T) at any T). 

The above argument is certainly correct for the non-interacting limit of the PAM, 
which is characterised strictly by simple one-electron states. Despite its superficial 
appeal however, it is incorrect beyond the confines of this limit: for it neglects completely 
[7] the effects of scattering arising from electron interactions, as embodied in the 
scattering rates r _1 (u;;T) = 7/(cj) (equation (3.3)) considered explicitly in §4.1; or 
equivalently (see equation (3.4)) in the imaginary part of the /-electron self-energy 
Tif(iv;T), which also controls the single-particle dynamics (§2,4). The effect of this 
many-body scattering on the optical conductivity may be inferred from equation (3.7a), 
from which it is straightforward to prove that if Ej(cu;0) is non-zero for |a>| > A, then 
Fhcl(w,0) is non-zero for uj > 2A. But the scale A above which £j(cj;0) effectively 
(HCL) becomes non-zero is of course the single-particle gap scale A g , whence (equation 
(6.2)) the optical gap corresponds to the indirect gap A ind . Four further points should 
be noted here, (a) We emphasize that this qualitative conclusion is not dependent on 
the LMA, although the latter indeed gives A ~ A g . In fact it arises also within IPT 
for which, as discussed in §5, £/((*;; 0) is non-vanishing for > 3A g (the central point 
again being that A oc A g ). (b) The above conclusion is certainly consistent (§5) with 
a transport gap (equation (5.1)) A tr ~ A g in the static conductivity: the latter is 
simply the uj = limit of the optical conductivity, and at the very least it is natural 
to expect the optical gap to be proportional to the low-temperature gap scale for static 
transport, (c) Recall that the scaling regime characteristic of the strong coupling Kondo 
lattice (where the quasiparticle weight Z — > 0) corresponds to any finite uj = uj/A g in 
the formal limit A g oc Z — > 0. But if the optical gap corresponded to the direct gap 
equation (6.3), this would require Ej(w;0) = for all < A dir oc Z 1 / 2 . Since such a 
frequency range clearly encompasses (and goes well beyond) the ^-scaling regime, this 
would entail vanishing scattering rates throughout the scaling regime; which is physically 
untenable, (d) Our discussion of the optical gap has in large part focussed implicitly 
on Fhcl(w,0) for the hypercubic lattice. But the same conclusion arises (essentially 
trivially) for the Bethe lattice Fbl(w, 0): from equation (3.7b) 

t 2 r° 

F bl {uj- 0) = - duu 1 D'iujD'iu! + uj) (6.4) 

^ J-UI 

where, since (§4) the conduction electron spectrum has a gap A g (measured from 
the Fermi level), the optical gap follows immediately as A ind = 2A g . This suggests 
in addition that we should expect the optical properties of the HCL and BL to be 
qualitatively similar, at least for sufficiently low frequencies and/or temperatures. That 
is indeed the case, as we show below. 

Figure 14 shows the resultant LMA optical conductivities F a (uj; T) (equations (3.7)) 
for both the BL and HCL, as a function of uj = uj/A g and for a range of temperatures 
T = T I A g . These are the universal forms, scaling in terms of uj and T with no 
dependence on bare material parameters [and by way of orientation we add that for 
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Figure 14. Optical conductivities in the scaling regime: Fbl(w;T) (left) and 
F hcl (uj- : T) (right) vs u/A g for f = T/A g =0 (solid), 0.5 (dotted), 1 (short dash), 
2 (point-dash), 5 (long dash) and 10 (double point-dash). The optical gap Aj„d = 2A„ 
is indicated in each case. 

gaps A g in the range ~ 10K — 300K, Co = 30 would correspond to frequencies uo in 
the range ~ 200 — 6000cm -1 ]. The indirect gap A in( i (= 2A g = 2ZV 2 /£*) is marked 
on the figures and is indeed seen to be the (T = 0) optical gap. The T = optical 
conductivity is 'sharper' in an obvious sense for the BL than the HCL, but in either 
case it is again the indirect gap that sets the scale for thermal destruction of the 
gap: by T = T/A g > 5 or so, the gap is well filled in and the optical conductivity 
essentially constant over the cu-range shown. These are of course characteristic features 
of experimental Kondo insulators [3-7], as discussed for specific materials in §7. Here we 
simply add that F a (oo; T)'s on the order of ~ 0.25 — 0.75 lead to absolute conductivities 
cr(u; T) = ^o~ F a (u; T) in the range ~ 2500 — 7500(ficm) _1 (taking a typical a (§3) of 
~ 3.10 4 (ficro) _1 ); values that are typical for Kondo insulators around room temperature 
[3-7]. 

Figure 14 also shows that on initially increasing T from 0, the optical conductivity 
acquires a Drude-like peak centred on uj = 0, which broadens with increasing T and 
is subsequently destroyed as the gap is filled in; as also seen in a QMC study of the 
HCL [15] (discussed further below). This is more clearly evident in figure 15, where the 
thermal evolution of the Drude peak for the HCL is shown for T = 0.4, 0.8 and 1.2. The 
figure also makes comparison to a Lorentzian form for the Drude peaks; the quality and 
range of which fit (and it is merely that) are self-evident. 

The optical conductivities shown in figure 14 refer to the universal cD-regime. To 
consider all w-scales, encompassing in particular frequencies on the order of the direct 
gap and beyond, bare model parameters must of course be specified. Figure 16 provides 
an example, for U/t* — 6 and V 2 jt\ = 0.2. For both the HCL and BL the optical 
conductivities F a (u;T) are shown for all uj over which F a {uj;T) is non-zero, as a 
function of Co = co/A g on a log scale from 10 _1 to 10 4 ; and for the same range of 
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Figure 15. Drudc-like peaks in the low-(<D,T) optical conductivity: Fhcl{^',T) vs 
co/Ag for T = T/ Ag—OA, 0.8 and 1.2. Lorentzian fits to the data (dotted lines) arc 
also shown. 



temperatures T = T/A g (up to 10) employed in figure 14. The insets in each case show 
the renormalized band structure cD±(e) = u±(e)/A g vs the free (V = 0) conduction 
band energies e/t*; with the branches w±(e) obtained from solution of (Re^(uj) =) 
Jr{u) = e. From this the fiducial direct gap A dir is obtained (occurring as expected 
for e = 0); it is indicated on the frequency axis, and for the particular chosen bare 
parameters is seen in either case to be ~ 10 2 times the indirect gap A ind (= 2A 9 ), or 
~ 0.05t*. Two features are immediately apparent in figure 16. First, unsurprisingly, 
that the optical conductivity is large on frequency scales on the order of the direct gap; 
particularly for the HCL where, as expected physically and known from previous work 
[15,21], Fhcl(u;T) is strongly peaked around u ~ A dir . 

The second point concerns the thermal evolution of the optical conductivity. The 
T = T/A g range shown in figure 16 corresponds to temperatures up to 5 times the 
indirect gap A ind = 2A 9 ; over which range, as shown in figure 14, the optical gap 'fills 
in'. As seen from figure 16 however, temperatures of this order have essentially no 
effect on the optical conductivity at frequency scales on the order of the direct gap 
A^r ~ 10 2 Aj nf j, which for all practical purposes retain their T = values. This is not 
surprising, for one should expect the optical conductivity on frequency scales uj ~ A^r 
to be thermally eroded only on temperature scales of the same order. This is indeed the 
case, as illustrated in figure 17 which (for the same bare parameters) shows the thermal 
evolution of F HC l(u;T) up to temperatures T = 200 (i.e. T ~ 0.7A dir ): significant 
thermal erosion on the direct gap scale sets by around T/A dir ~ 0.1 — 0.3, and is well 
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Figure 16. Optical conductivities on all w-scalcs: Fhcl(^' : T) (top panel) and 
Fbl{w, T) (lower panel) vs uj/A 9 on a log-scale, for T = T/A g —0 (solid), 0.5 (dotted), 
1 (short dash), 2 (point-dash), 5 (long dash) and 10 (double point-dash). The 
indirect and direct gap scales are also indicated. Temperatures in the range shown 
produce essentially no effect at w's on the order of the direct gap scale A-dir- Insets: 
renormalized band structure u±(e)/A g vs the free (V — 0) conduction band energies 
e/t*. 

established at the highest temperature shown in the figure. 

It is clear from the above that the indirect and direct gap frequency scales, each 
determined by but with their different dependences upon the quasiparticle weight Z 
(equations (6.2,3)), are qualitatively distinct in the strong coupling Kondo lattice regime 
where Z <C 1; as too in consequence are the corresponding temperatures for which 
dynamics on these respective scales evolve. But this 'clean' separation of scales will not 




Figure 17. For the same strong coupling parameters as figure 16, Fhcl{oj\ T) vsui/ A g 
for a much wider range of temperatures T — T/A g up to 200 (i.e. T I A dir ~ 0{1)); 
significant thermal erosion on the direct gap w-scale occurs only for temperatures of 
the same order. 



of course be evident in the (ou, Independences of the optical conductivity if one is e.g. 
restricted technique- wise either to low interaction strengths (where Z ~ 0(1)) or to high 
temperatures. This is the case with Quantum Monte Carlo, and is we believe the reason 
why the above scale separation is not apparent in the QMC work of [15]. This lack 
of scale separation is also intrinsic to IPT [21], since its resultant quasiparticle weight 
decays algebraically rather than exponentially in the interaction strength. Numerical 
renormalization group calculations [12,13] by contrast are not constrained in the above 
sense, although NRG calculations of finite-T dynamics and transport properties of the 
PAM have not to our knowledge been performed thus far. 



7. Experiment 

Experimentally, Kondo insulators have been widely studied via an impressive range 
of techniques (for reviews see e.g. [3-7]). Here we consider briefly three prototypical 
materials among those for which the most extensive and reliable data is available, viz 
Ce 3 Bi4Pt 3 , SmB 6 and YbB 12 . Each naturally possesses features specific to itself, but 
the well known commonality of behaviour between the different materials is of course 
the dominant theme [3-7]. Our aim here is simple: to compare the present theory 
directly to experimental results for d.c. transport and optical conductivities. We also 
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emphasise that subsequent comparison to experiment does not involve multi-parameter 
fits, or a detailed knowledge of the bare material parameters; and in these terms is 
minimalist. The essential point of scaling, as detailed in the preceding sections, is 
that in strong coupling the temperature dependence of transport/optics is encoded in 
T = T/Ag alone, independently of the bare model parameters. The gap A g (oc A tr ) 
may thus itself be obtained from experiment (as below); given which, the theory then 
predicts the full T-dependence. The same comment applies to the (^-dependence of the 
optical conductivity, which in the scaling regime depends only on to = oj/A g (albeit 
of course that 'non-universal' cj's (§6) are also accessed experimentally). All this is 
naturally based on the assumption that the experimental materials are indeed strongly 
correlated, which we take for granted unless there appear to be experimental hints to 
the contrary (as we suggest below may be the case for YbB^). 

Granted the above there is one unknown in the theory, viz the constant o"o entering 
the conductivity a(cu;T) = |cr F a (u;; T) as an overall scale factor (see e.g. equation 
(3.6)); and which, as noted in §3, should realistically lie in the range ~ 10 4 — 10 5 (ficm) _1 . 
In comparing to experiment we take a as a parameter, which in practice is indeed found 
to lie in the above range; although even this is not strictly necessary (knowledge of <To 
may be bypassed entirely by taking e.g. the room temperature d.c. conductivity as a 
reference and considering ct(uj; T)/<x(0; 300K)). One must of course choose whether to 
compare experiment to theoretical results arising from the HCL or BL, the two host 
lattices we have considered explicitly. Here we simply make comparison to whichever 
of the two appears optimal, which in practice means the hypercubic lattice for all but 
Ce 3 Bi i Pt 3 ; a more realistic free (V = 0) conduction band density of states Po( e ) could 
readily be employed, but the comparisons below suggest this would be barely necessary. 

7.1. Ce 3 Bi 4 Pt 3 . 

Experimental results for the resistivity p(T) of this classic Kondo insulator are shown in 
figure 18. These are compared to results arising for p(T) (= l/er(0; T)) from the present 
theory for the BL, taking a ~ 4.2 x 10 4 (ficm) _1 . The transport gap (equation (5.1)) 
inferred experimentally is A tr ^ 35K [48], from which (§5) A g = A tr follows. 

The quantitative agreement between experiment and theory is clear from figure 18, 
being excellent for T > 50K or so, as seen in particular from the d.c. conductivity shown 
directly in the inset. The agreement extends up to T = 300K or T = T/A g ~ 10, 
an appreciable multiple of the gap where (see e.g. figure 10) the high-temperature 
logarithmic asymptotics of the Kondo regime are being approached; so that the 
temperature range shown spans essentially the full range of expected physical behaviour. 
The activated insulating nature of the low-T transport (equation (5.1)) is likewise 
evident in the figure, although transport in CesBi^Pts at the lowest temperatures is 
dominated by variable range hopping between extrinsic states in the gap [48] which the 
theory does not of course seek to include. 

As discussed in §6, the d.c. transport and optics should be consistent in the sense 
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Figure 18. Ce^BiiPt^ resistivity p{T) (in Qcm) vs T up to 300K. Experimental 
results [48], open circles; theory, solid line. Inset: corresponding results for the d.c. 
conductivity <r(0;T). 

that the (T = 0) gap in the optical conductivity should correspond to the indirect 
gap, A ind = 2A g ; and the thermal evolution of the optical conductivity should likewise 
be controlled by the indirect gap scale. Taking ~ 70K ~ 50cm -1 from the d.c. 
experiments as above, figure 19 (top panel) shows the theoretical optical conductivity 
a{uj;T) = ^o~qFbl{w,T) vs uj up to 1000cm -1 , for temperatures T = 0,25,50,75,100 
and 300-fT (each curve thus corresponding to a particular 'realisation' of the universal 
BL optical conductivities shown in figure 10). The lower panel in figure 19 shows 
corresponding experimental results [49] for uj > 50cm -1 , obtained from Kramers-Kronig 
transformation of reflectivity measurements; we have also indicated the indirect gap 
on the a>axis, and have marked experimental values of the d.c. conductivity [48] on 
the vertical axis. The level of agreement between experiment/theory for the optics is 
self-evident from figure 19, and we regard it as rather good. The relevant temperature 
scale indeed appears to be the indirect gap [7,50]: as noted in [49], for T between 
100K ~ lAA in( i and 300/T the optical gap is well filled in and the optical conductivity 
nearly constant in the far infrared; and the gap begins to develop markedly only below 
T ~ 75K ~ A in d (which we add is naturally consistent with thermal evolution of the 
single-particle dynamics, see e.g. figure 4). 

As seen from the experimental data in figure 19, the low-temperature d.c. 
conductivities are larger than their a.c. counterparts at the lowest uj = 50cm _1 , the 
difference diminishing with increasing T and being barely perceptible by T ~ 75K ~ 



Dynamics and transport properties of Kondo Insulators 
— i 1 1 1 1 1 1 1 1 — 



40 




Figure 19. Ce^Bi^Ptz optical conductivity cr(w;T) (in Q,' 1 cm' 1 ) vs lo up to 1000 
cm' 1 , and for temperatures T=25,50,75,100 and 300K (in obvious sequence). Top 
panel: theory (including T = 0); bottom panel: experimental results for uj > 50cm' 1 
[49], with experimental d.c. conductivities marked by crosses on the vertical axis. The 
optical/indirect gap Ai nc i is indicated; for 'A c ', see text. 



A in d- This we believe to be symptomatic of the Iow-cj Drude-like peak discussed in 
§6 (figures 14,15) that lies below the detection limit of 50cm _1 but is apparent (albeit 
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weakly) in the theoretical o~(ou; T); and further evidence for which arises in the systems 
discussed below (see figures 21,23). Finally, a remark on the optical gap itself. As noted 
in [49], linear extrapolation to zero of the steep part of the experimental cr(u; T = 25K) 
would suggest a gap value on the order of 100K. This is certainly consistent with 
Ai n d — 70K ~ 50cm _1 , inferred as above from d.c. transport. On the other hand a 
charge gap A c ~ 300cm" 1 has been identified in [49], for natural reasons evident in the 
experimental data shown in figure 19 (and with A c indicated on both panels in figure 
19). Our point here is simply that A in d and A c are fundamentally equivalent scales (each 
being proportional to the quasiparticle weight Z). The former is 'correct' insofar as the 
optical gap is strictly a T = notion; while the latter is natural if one wishes instead to 
focus on the incipient development of a gap coming from the 'high'-temperature regime 
T > A ind . 

7.2. SmB 6 . 

Samarium hexaboride provides another long studied [51], prime example of a Kondo 
insulator. Here we refer to a recent comprehensive study [52] of both static transport 
and low-energy electrodynamics, performed on a high quality single crystal sample; 
with the optical conductivity obtained directly via sub-millimeter spectroscopy in the 
frequency range 5 — 36cm -1 [52], and by Kramers- Kronig analysis of reflectivity in the 
infrared [52,53]. We restrict our considerations mainly to temperatures T > 8K, since 
for lower temperatures variable range hopping again arises. 

Experimental results for the resistivity p(T) [52], spanning five orders of magnitude, 
are shown in figure 20; and compared to theoretical results for the hypercubic lattice 
(as one might anticipate to be appropriate for a clean sample), taking cr = 4.7 x 
lO A (Qcm)~ 1 and A g = 101K. Between T = 8K and 25K the experimental d.c. 
conductivity/resistivity has the activated form equation (5.1), with a transport gap 
A tr ~ 3.5meV (~ 41i^) [52]. As discussed in §5, for the hypercubic lattice we find 
A tr = bA g with the constant b = 0.40; hence A g ~ 8.75meV ~ 101/T as above. The 
main figure shows p(T) on a log-scale vs T; while the insets show the corresponding 
d.c. conductivity vs both 1/T (to exemplify in particular the activated regime) and 
temperature on a log-scale. The agreement between theory/experiment for T > 8K is 
rather striking; for T < 8K the experimental variable range hopping behaviour is of 
course evident from the right inset to the figure. 

The consequent optical conductivity c(c<j; T) for the HCL is shown vs u in figure 
21 (top panel) on a log-log scale, for temperatures T = 3, 13, 16, 18 and 300JT. 
Corresponding experimental results for the same temperatures [52] are shown in the 
lower panel, including the extrapolated conductivities (dashed lines) obtained from 
the phenomenological fit to the data employed in [52]. The cu-range shown, up to 
~ 3 x 10 4 cm _1 , naturally encompasses non-universal scales at high-frequencies, and the 
calculations were performed specifically for U/t* — 4.5 and V 2 /tl = 0.2. The system 
is however strongly correlated for these parameters, so we emphasise that the resultant 
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Figure 20. SttiBq resistivity p(T) (in ficm) vs T up to 300K. Experimental results 
[52], dashed line; theory, solid line. Left inset: corresponding d.c. conductivity a(0:T) 
vs T on a log-log scale. Right inset: er(0; T) on a log scale vs 1/T; for T < 8K variable 
range hopping arises experimentally. 

'low '-a; conductivity (up to ~ 2000cm -1 or so in practice) lies in the u = u/A g scaling 
regime that is actually independent of the bare parameters (§'s 4-6): the choice of 
bare parameters 'matters' only at the high frequency end, and we have simply chosen 
those above as illustrative. We would however add that for the chosen parameters the 
quasiparticle weig ht Z ~ 5xl0~ 3 , which is in qualitative accord with the experimentally 
inferred effective mass (1/Z ~) m*/m ~ 10 2 (from [52]) and ~ 10 3 from a previous 
study [54] (and which values themselves attest to the correlated nature of SmB 6 ). 

The first point to note here is that the theoretical optical gap A ind = 2A g , indicated 
on figure 21, is A ind ~ 17.5meV ~ 200K (from A g obtained as above). This accords 
remarkably well with the gap of 19 ± 2meV inferred experimentally from the optical 
conductivity [52], and we emphasise that there is no a priori connection between these 
two ways of obtaining the optical gap - - as above, the theoretical optical gap is 
inferred directly from the much smaller transport gap A tr . There is in consequence no 
conundrum between a transport gap of ~ 40i^ and an optical gap of ~ 200i^. The level 
of agreement between theory/experiment is self-evident in figure 21, and encouraging 
both in terms of its cu-dependence and thermal evolution. One small twist may also 
be added. For the phenomenological fitting employed in [52], it was found that for a 
complete description of the conductivity spectra an additional parameter (termed l (T mj „') 
had to be introduced in the form of an additive, frequency-independent contribution to 
the optical conductivity; with a min ~ 12 {Vtcm)^ 1 in the low-temperature regime [52]. 
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Figure 21. SmB 6 optical conductivity a(w;T) (in fl~ 1 cm~ 1 ) vs w on a log-log 
scale, for temperatures T — 3,13,16,18 and 300K. Top panel: theory (with the 
theoretical optical/indirect gap indicated). Inset: corresponding results when a 
constant 12Q _1 cm _1 is added. Bottom panel: experimental results [52]. Circles are 
from sub-mm data [52], solid lines from reflectivity spectrum via Kramers-Kronig (KK) 
[53]. Error bars refer to IR conductivity obtained from KK analysis of the T = 3K 
reflectivity spectrum assuming 0.5% uncertainty [52]. The shaded area corresponds to 
the experimental optical gap [52] . The short dash lines show the extrapolated fit used 
in [52], and the arrows indicate the experimental d.c. conductivity. 



Without wishing to speculate here on the origin of the a m i n , the inset to figure 21 
(top panel) shows the effect of simply adding a constant 12 (Hem)" 1 contribution to 
our theoretical optical conductivity. For low-temperatures its effect is noticeable in the 
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~ 10 — 100cm" 1 range, particularly at the lowest T = 3K; and the improved agreement 
with experiment is clear. 

Finally we mention that we have taken no consideration above of the presence, in 
the Kondo insulating gap, of an additional narrow donor-type band which is known to 
occur experimentally [52] in SttiBq (and may possibly be due to impurities). While 
certainly of interest in itself, the present theory has of course nothing to say about it; 
and it plays no role in the extent to which, as above, theory concurs with experiment. 



7.3. YbB 12 . 

Resistivity [55] and optical measurements [56] have likewise been performed on single 
crystal YbBu, the only known Fft-based Kondo insulator. Experimental results for 
p(T) up to ~ 350/T are shown in figure 22, and for 15K < T < AOK exhibit activated 
behaviour (equation (5.1)) with a transport gap A tr ~ 68K. Optical conductivity 
results [56], again obtained via Kramers- Kronig from reflectivity spectra combined with 
a Hagen-Rubens extrapolation at low-u;, are shown in the top panel of figure 23. The 
experimental optical gap is determined as Aj„d ~ 25meV = 290K [56] ; while the strong 
IR peak around ~ 0.2 — 0.25eV" as naturally interpretable [7] in terms of direct gap 
excitations (see e.g. figures 16,17 for the HCL). 
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Figure 22. YbB\2 resistivity p(T) (in Hem) vs T. Experimental results [55], open 
circles; theory, solid line (for parameters specified in text). Inset: corresponding d.c. 
conductivities vs T. 



It is the temperature dependence of the experimental optical data which suggests to 
us that YbB\2 may be in an intermediate-weak coupling regime. As seen from figure 23 
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Figure 23. YbB 12 optical conductivity ct(cj;T) (in Q~ 1 cm~ 1 ) vs w (in eV), for 
T=20,78,160 and 290K (in obvious sequence). Top panel: experimental results [56]. 
Bottom panel: theory (for same parameters as figure 22). 

[56], increasing temperatures up to T = 290K = Ai n d leads to significant redistribution 
of spectral weight at much higher energy scales on the order of the direct gap and beyond. 
This is not behaviour characteristic of strong coupling, as evident from the discussion 
of §6 (figures 16,17). It is however typical of intermediate- weak coupling interactions, 
theoretical consideration of which thus requires specification of bare model parameters. 
In the following we consider specifically U/t* = 1.65 and V 2 jt\ = 0.2: these values 
should not of course be taken too seriously per se, but they lead within our approach 
to behaviour representative of intermediate-weak coupling and should be viewed simply 
as such. 

We find with the latter that the low-temperature transport has as expected the 
activated form equation (5.1), with A tr = bA g and the constant b = 0.47 (as opposed 
to b = 0.40 for the HCL in strong coupling). Equating the theoretical A tr with the 
experimental A tr ~ 68K, and taking a ~ 10 5 (f2cm) -1 , direct comparison between 
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the T-dependence of the theoretical and experimental p(T)'s is shown in figure 22; 
and the agreement with experiment is seen to be rather good across essentially the full 
temperature range. Does this lead to a consistent description of the optical conductivity? 
The answer is yes in two senses. First, the theoretical estimate of the optical gap 
Ai n rf = 2A g = (2/6) A tr follows directly as A in( i ~ 290K = 25meV using only the 
experimental transport gap above; which value coincides with A ind obtained directly 
from the optical experiments (a situation analogous to that for SmB 6 discussed above) . 
Second, the resultant theoretical optical conductivity is shown in the bottom panel of 
figure 23. It is seen to accord well qualitatively with the experimental results (top 
panel), both in terms of its overall cj-dependence and thermal evolution; including the 
redistribution of spectral weight on energy scales beyond the direct gap, for temperatures 
up to T = 290K that merely corresponds to the indirect gap itself. (Improved agreement 
with experiment could no doubt be obtained by playing with the bare parameters, but 
would add little new.) 

8. Summary. 

We have developed in this paper a non-perturbative local moment approach to 
dynamics and transport properties of the symmetric periodic Anderson model, the basic 
microscopic model for understanding small-gap Kondo insulator materials [3-7]. Our 
primary focus has naturally been the strong coupling, Kondo lattice regime. Here the 
system is characterised by the low-energy scale A g which, being exponentially small in 
strong coupling, leads to a pristine separation between low- and high-energy scales; and 
hence to 'universal scaling' of dynamics/transport in terms of Co = uo / A g and T = T/A g 
alone, with no explicit dependence on bare model parameters. It is this single, indirect 
gap scale A g that is of paramount importance in controlling the physical properties of 
the system that we have investigated systematically here; for it determines the single- 
particle spectral gap, the transport gap for d.c. conductivity and the optical gap in 
the dynamical conductivity, all of which are simply proportional to each other. It sets 
the scale for thermal evolution of single-particle dynamics and d.c. transport, from the 
gapped/activated behaviour symptomatic of the low-temperature insulator through to 
the incoherent single-impurity physics that is found to arise naturally for T ^> 1. And 
likewise it is A g that sets the thermal scale for 'filling' the optical gap with increasing 
temperature; the much higher direct gap scale naturally also being apparent in the co- 
dependence of the optical conductivity, but in strong coupling being essentially irrelevant 
as a thermal scale for its evolution. 

Notwithstanding the innate simplicity of the PAM itself, and the range of material- 
specific factors it naturally omits, the present theory also appears to account well 
for experiments on materials such as C e^Bi^Pt^,, SmB 6 and YbB 12 ; with many 
characteristic features arising theoretically apparent in experiment, and a mutually 
consistent picture of d.c. transport and optics arising. We believe this lends further 
support to the essential veracity of both the underlying model (the PAM within DMFT) 
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as well as the present theory, the local moment approach. Further development of the 
LMA to encompass the asymmetric PAM, and hence heavy fermion metals, will be 
reported in subsequent work. 
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Appendix. 

Here we sketch the steps leading to equation (5.8) for the leading high-T behaviour 
of p BL (T) for the Bethe lattice. Using equation (3.7b) for uo = 0, together with 
< D c (e\uo) > t = D c (uj) (the c-electron spectrum), F BL (Q;T) is given generally by 

Fbl(0; T) = tl jH du - d -^l [D c (u)\ 2 . (1) 

For the BL, the Feenberg self-energy S(uo) = \tlG c (uo)\ so equation (2.5a) implies 

G» = [ 7 (o;) - \tlG\u)]- 1 (2) 

which determines the 7-dependence of G c = G c [7]. For T ^> 1 in the scaling regime, 
I7I <C £*; which corresponds physically to scattering rates (7/(0;) = 7rp 7/(^) =) 
f _1 (u;;T) <C 1 for (all) finite |u)|. Using equation (2), the leading asymptotics of 
D c {u) = — ^ImG c {ui) is given for | 7 | <C 1 by 

D c (u) ~ p [l - f P07/H + = po[l - T) + 0(f- 2 )] (3) 

whence from equation (1): 

F BL (0-T) ~ [p uf{\ + J°° du r~\u;T)} (4) 

Equations (5.6,7) for f _1 (c<j;T) then lead to 



where f(y) = [e y + l]" 1 ; and using L(y;f) -> 1 as f -> 00 (§5) gives F BL {0;T) ~ 
[po^] 2 {l- 16 fa 2 2 (:f) }- Using this in equation (5.4) for p' BL (T) (together with [p t*] 2 = ^), 
gives directly the leading large-T behaviour 

PBLl J ~ 16/n 2 (f) 1 J 

as sought. 
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